UNCLASSIFIED 


..  401  770 

R«fVtoduc*d 

if  tlu 

ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


oo 


FINAL  RIPORT  NO.  63-94 
This  report  covers  the  period 
August  1, 1961  to  July  31. 1962 


sO/> 

% 


% 


STUDY  OF  ELECTRICAL  BREAKDOWN  CONDITIONS 
IN  THE 

AERODYNAMIC  FLOW  FIELD  OF  A  HYPERSONIC  VEHICLE 


by 

Morton  Rudin 
Boris  Regent 
Charles  E.  Noble,  Jr. 


prepared  for 

OFFICE  OF  NAVAL  RESEARCH.  WASHINGTON.  D.C. 
and 

ADVANCED  RESEARCH  PROJECTS  AGENCY.  WASHINGTON.  D.C. 


cb 

qc 3 


<c  ~> 

c_3  <c 


RESEARCH 


Contract  No.  Nonr  3581(00) 
Under  ARPA  Order  No.  207-62 
Project  Code  No.  7400 
"Microwaves  in  Hot  Gases" 

January  30,  1963 
Vidya  Project  No.  64/C 


DEVELOPMENT 


■.  1)  1  i  A 


i if  i-5 » 


j  r^j  -l  J  jjJj  j 


TISIA 


A  DIVISION  or 


CORPORATION 


VIDYA 


7  he  uord  Vidyu,  t  it  ken  from  the  Vedanta  philosophy  of  the  Hindus, 
means  knowledge.  The  symbol  used  to  denote  the  I  idya  organization 
is  the  li  tter  "Y"  from  Sanskrit,  the  ancient  language  of  India. 


MECHANICS  ...  .  PHYSICS...  .ANALYSES 


VIDYA  COPY  NO. 


PINAL  REPORT  NO.  63-94 
Thli  report  covet*  tho  period 
Aufuit  1,  1981  to  July  31,  1962 

STUDY  OF  ELECTRICAL  BREAKDOWN  CONDITIONS 
IN  THE 

AERODYNAMIC  FLOW  FIELD  OF  A  HYPERSONIC  VEHICLE 

Morton  Rudln 
Borii  Rogont 
Chariot  E.  Noble,  Jr. 

prepared  for 

OFFICE  OF  NAVAL  RESEARCH,  WASHINGTON,  D.  C. 
end 

ADVANCED  RESEARCH  PROJECTS  AGENCY,  WASHINGTON,  D.  C. 


Contract  No.  Nonr  3581(00) 
Under  ARPA  Order  No.  207-62 
Project  Code  No.  7400 
"Microwave*  In  Hot  Qetei” 


January  30,  1963 
Vitlya  Project  No.  64/C 


Thl*  teteerch  ie  e  pert  of  Project  DEFENDER  tpontored  by  the  Advanced 
Reoaotch  Project*  Agency,  Department  of  Defen**,  end  technically  edmlnletered 
by  the  Fluid  Dynamic*  Branch  of  the  Office  of  Nival  Retearch. 


Reproduction  In  whole  or  In  pert  I*  permitted  for  any  purpott  of 
The  United  State*  Government 


V  I  D  Y  A 


A  DIVISION  OF 


CORPORATION 


1450  PAGE  MILL  ROAD  #  PALO  ALTO,  CALIFORNIA 
TILt  DAVINPORT  1-2455  TWXs  415  492-9270 


ABSTRACT 


A  numerical  technique  has  been  formulated  for  solving  the 
antenna  breakdown  equation  in  the  presence  of  gaseous  flow  for 
typical  antenna  electric  field  and  aerodynamic  flow-field  con¬ 
figurations.  The  resultS'Of  computations  for  nonflow  conditions 
agree  well  with  existing  data.  Breakdown  conditions  were  computed 
for  sample  cases  of  antenna  type,  vehicle  shape,  and  velocity- 
altitude  combinations.  The  results  in  air  show  that  once  anibi- 
polar  diffusion  has  been  established  increasing  the  vehicle  velocity 
at  a  fixed  altitude  usually  increases  the  breakdown  field  strength. 
This  increase  is  found  to  be  caused  principally  by  the  increase  of 
gas  density  near  the  antenna  and,  to  a  much  lesser  degree,  by  the 
effect  of  the  boundary-layer  flow  velocity  in  sweeping  electrons 
out  of  the  high  electric  field  regions . 

The  kinetic  theory  concepts  involved  in  the  above  calculations 
have  been  critically  examined  and  the  gaseous  material  properties 
necessary  for  these  calculations  have  been  estimated.  The  inviscid 
and  viscous  flow  fields  in  which  the  test  antenna  is  immersed  were 
calculated  for  wedges  and  cones,  using  equilibrium  real -gas  theory, 
and  the  results  of  these  calculations  were  used  as  inputs  for  the 
previously  discussed  calculations. 

Finally,  the  results  of  these  calculations  are  summarized  and 
conclusions  and  recommendations  for  future  work  are  given. 
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STUDY  OP  ELECTRICAL  BREAKDOWN  CONDITIONS  IN  THE 
AERODYNAMIC  FLOW  FIELD  OF  A  HYPERSONIC  VEHICLE 

1 .  INTRODUCTION 

This  report  covers  work  done  on  Contract  Nonr  3581 (00) , 

Project  No.  061124,  during  the  period  1  August  1961  to  1  August 
1962.  The  purpose  of  this  project  is  to  perform  a  "Study  of 
Electrical  Breakdown  in  the  Hypersonic  Flow  Field  of  a  Hypersonic 
Vehicle . " 

The  difficulties  of  propagating  electromagnetic  energy  through 
highly  ionized  gases ,  such  as  may  occur  in  a  plasma  sheath  sur¬ 
rounding  hypersonic  vehicles,  have  been  much  discussed.  The  atten¬ 
uation  of  weak  signals  in  such  plasmas  is  fairly  well  understood 
theoretically,  and  for  such  a  situation  the  electron  random-energy 
distribution  or  random-velocity  distribution  is  very  little  dis¬ 
turbed  by  the  presence  of  the  electromagnetic  field,  and  calcula¬ 
tions  based  upon  an  unperturbed  distribution  function  usually 
suffice.  Contrasted  with  this  case  is  that  when  a  strong  electro¬ 
magnetic  field  is  present  in  the  gas,  such  as  occurs  near  an  antenna 
surface.  In  this  case  the  presence  of  the  field  strongly  perturbs 
the  initial  electron  random-velocity  distribution  function  by  adding 
energy  to  the  electrons  and  may,  under  certain  conditions,  lead  to 
very  high  ionization  rates,  a  condition  known  as  "breakdown."  This 
high  ionization  surrounding  an  antenna  may  change  the  radiation 
characteristics  of  the  antenna  and  drastically  impair  propagation 
of  energy  from  the  antenna. 

For  electrical  breakdown  in  the  flow  field  of  a  hypersonic 
vehicle,  we  are  interested  in  the  breakdown  condition  for  both 
pulsed  and  c-w  electromagnetic  waves  in  hot  air  (possibly  not  in 
thermodynamic  equilibrium  and  with  ablation  additives)  and  with 
the  gas  flowing  in  an  open  geometry  over  a  surface  containing  the 
source  of  the  electromagnetic  field.  The  flow  field  is  not  uniform 
and  the  alternating  electric  field  is  not  of  uniform  amplitude. 

All  of  these  conditions,  as  examination  in  this  report  shows,  can 
and  do  have  some  effect  upon  the  breakdown  conditions. 
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Experimental  studies  of  the  breakdown  condition  in  the  past 
include  c-w  breakdown  of  a  gas  in  a  closed  container  (data  and 
experiments  reviewed  in  Refs.  1,  2,  and  5),  pulsed  a-c  breakdown 
in  a  closed-container  (Ref.  3),  no-flow  pulsed  a-c  breakdown  in 
an  open  region  over  a  surface  containing  the  antenna  (Ref.  4) , 
and  breakdown  with  flow  between  plates  or  within  tubes  (Ref.  6) . 

All  of  these  were  for  relatively  cold  gases  (approximately  300°  K) . 
In  addition,  there  have  been  measurements  of  breakdown  field 
strengths  in  flames  of  uncertain  composition  (Refs.  7a  and  7b)  and 
in  a  plasma-filled  cold-air  region  (Ref.  7c) . 

Theoretical  breakdown  studies  in  the  past  include  those  for 
closed  containers  with  c-w  (Refs.  1,  3,  and  5)  or  pulsed  (Ref.  3) 
uniform  or  nonuniform  electric  fields.  There  is  also  a  closed 
solution  for  breakdown  in  a  gas  flowing  with  constant  velocity  in 
a  tube  or  between  plates  with  an  applied  electric  field  constant 
in  a  finite  region  and  zero  elsewhere  (Ref.  8) . 

All  of  these  studies  were  not  concerned  with  the  effects  of 
hot  (rather  than  cold)  air.  Also,  of  all  the  papers  cited  above, 
only  one  (Ref.  4)  has  a  surface  geometry  and  electric  field  geom¬ 
etry  pertinent  to  the  problem  with  a  reentry  vehicle.  In  addition 
to  all  the  papers  mentioned  above,  there  are  some  cursory  attempts 
to  compute  breakdown  conditions  for  reentry  problems  (Refs.  9,  10, 
and  11).  In  Reference  9,  the  breakdown  field  strength  is  estimated 
for  air  (not  flowing)  in  a  field-containing  region  which  is  assumed 
to  be  ~K/2  long.  The  problem  is  first  treated  as  diffusion- 
controlled  breakdown  with  diffusion  length  V 2.  Subsequently, 
the  problem  is  treated  as  ionization  growth  in  a  uniform  field  for 
finite  time,  without  diffusion  and  with  all  electrons  produced  in 
previous  pulses  swept  out.  In  Reference  10  these  estimates  are 
modified  to  include  ambipolar  phenomena  and  "initial  ionization 
level . " 

We  find  (as  will  be  shown)  that  this  last  analysis  is  not 
applicable  to  our  problem  for  the  following  reasons  -  in  addition 
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to  the  unqualified  use  of  ionization  rates  for  cold  air.  First, 
we  find  that  diffusion  into  the  antenna  containing  wall  contributes 
appreciably  and  may  even  predominate  over  the  diffusion  across  the 
field.  Second,  we  find  that  no  initial  ionization-level  correction 
is  necessary  for  the  multiple  pulse,  diffusion-pulse-width  con¬ 
trolled  case  and  that  such  a  correction  only  has  meaning  when  one 
is  trying  to  inhibit  breakdown  by  use  of  short  pulse  width  alone. 
Third,  it  is  difficult  to  be  sure  that  the  electrons  from  one  pulse 
are  completely  swept  out  before  the  advent  of  the  next  pulse,  since 
the  boundary-layer  velocity  is  zero  at  the  wall.  One  must  treat 
the  problem  in  a  completely  coupled  manner. 

In  Reference  11 ,  the  reentry  problem  is  again  treated  as  ioni¬ 
zation  growth  in  a  uniform  field  during  a  finite  time.  The  only 
major  modification  over  References  9  and  10  is  that  the  time  con¬ 
cerned  is  taken  as  the  smaller  of  V2V  and  T,  where  V2V  is 
the  time  to  traverse  the  assumed  width  of  the  field  when  flow  is 
at  velocity  V,  and  T  is  the  pulse  width.  Diffusion  and  multiple 
pulses  are  not  considered.  Also,  the  data  used  are  again  for  cold 
air. 

We  treat  the  problem  in  a  completely  coupled  manner  and  formulate 
a  procedure  by  which  any  velocity  variation  and  any  nonuniform  field 
may  be  considered  in  the  breakdown  computations  for  the  desired 
geometry.  These  concepts  are  discussed  generally  in  Section  2, 
but  are  not  fully  formulated  for  the  problem  at  hand  until  Section  5. 

In  Section  3,  restrictions  on  the  applicability  of  experimental 
data  obtained  in  various  laboratory  situations  are  analyzed  and 
computed.  The  point  is  that  such  data  are  not  always  applicable 
to  pulsed  or  flowing  problems,  even  in  cold  air.  Since  the  micro¬ 
scopic  processes  governing  ionization  growth  and  breakdown  phenomena 
are  dependent  on  the  velocities  of  the  electrons  involved  (by  col¬ 
lision)  in  these  processes,  it  became  necessary  to  study  the  random- 
velocity  distribution  functions  for  the  electrons  and  their  rates 
of  relaxation  under  the  physical  conditions  imposed  by  the 
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environment.  In  particular,  the  macroscopic  quantities  which  must 
enter  into  the  evaluation  of  breakdown  conditions ,  such  as  electron 
production  and  loss  rates ,  are  those  values  for  these  quantities 
which  have  been  averaged  over  the  pertinent  distribution  function. 
All  the  available  data  on  these  macroscopic  quantities  are  from 
measurements  made  under  conditions  which  involved  one  particular 
type  of  electron  random-velocity  distribution  function;  therefore, 
caution  must  be  exercised  in  utilizing  these  data  under  physical 
conditions  which  are  vastly  different  from  those  under  which  they 
were  determined  or  in  which  transient  effects  are  present. 

In  Section  4,  the  velocity  profiles  in  reentry  boundary  layers 
are  determined  for  two  sample  body  shapes  and  typical  flight  condi¬ 
tions  and  matched  with  a  function  for  computation.  Other  pertinent 
boundary-layer  properties  are  also  deduced. 

To  test  the  breakdown  computation  procedure,  computations  are 
made  in  Section  6  at  zero  velocity  to  compare  with  an  analogous 
experiment  in  the  proper  geometry.  In  Section  7,  computations  with 
a  velocity  profile  similar  to  the  reentry  boundary  layer  for  cones 
and  wedges  are  conducted.  Section  8  contains  a  summary  of  the  work 
accomplished  in  this  report,  conclusions  reached,  and  recommenda¬ 
tions  for  future  work. 
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2.  THE  ELECTRON  CONSERVATION  EQUATION  AND  THE  CONCEPT  OF  BREAKDOWN 

2.1  The  Electron  Conservation  Equation 

The  purpose  of  this  study  Is  the  determination  of  the  electric 
field  strength  necessary  to  cause  breakdown  In  a  gas  steadily  flow¬ 
ing  along  a  surface  and  through  an  applied  electric  field.  The 
flow  velocity  is  assumed  to  be  in  a  uniform  direction  parallel  to 
the  surface  and  its  magnitude  to  vary  with  distance  away  from  the 
surface.  The  electric  field  is  also  a  function  of  the  coordinates 
and  should  resemble  the  pattern  of  an  antenna  positioned  within 
the  surface.  Take  coordinates  (see  Fig.  2.1)  such  that  the  flow 
velocity  V  is  in  the  direction  of  the  z-axis,  the  surface  S  is 
in  the  x-z  plane,  and  y  is  distance  above  the  surface  S. 

The  equation  governing  the  electron  density  n(r,t)  is  then 

_  V  •  +  V2  (Dn)  +  (vj  -  v A  -  vR)  n  (2.1) 

where  D  is  the  diffusion  coefficient,  V  is  the  flow  velocity, 
and  (vj  -  vA  -  vR)  is  the  net  ionization  frequency  defined  in 
the  list  of  terminology.  The  origin  of  this  equation  is  discussed 
in  References  12,  13,  and  1;  it  is  also  discussed  in  Section  3 
with  Appendix  A.  Discussed  immediately  below  are  some  of  the  more 
pertinent  assumptions  involved  in  the  form  of  this  equation. 

2.1.1  The  linear  electron  production  term  and  boundary  conditions 

A  common  argument  for  the  linearity  of  the  electron  production 

term,  Vjn,  where  Vj  is  the  ionization  frequency  averaged  over 
the  velocity  distribution  function,  may  be  made  as  follows: 

The  electron  production  is  proportional  to  n  only  if  ioni¬ 
zation  produced  by  collisions  among  the  heavier  ions  and  neutrals 
is  negligible  compared  to  ionization  produced  by  electron  colli¬ 
sions  with  heavier  particles.  Since  we  will  be  interested  in 
determination  of  the  "breakdown  conditions,"  since  the  probability 
of  ionization  by  electron  impact  is  several  orders  of  magnitude 
greater  than  by  heavy  particle  impact,  and  since,  in  addition,  the 
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electrons  will  be  more  energetic  than  the  heavier  particles  once 
they  are  in  the  applied  field,  the  ionization  induced  by  heavy 
particle  impacts  may  usually  be  neglected  so  long  as  some  elec¬ 
trons  are  on  hand  initially. 

The  above  argument  is  unnecessary  since  the  ionization  due 
to  collisions  among  the  heavy  particles  is  actually  not  neglected 
by  use  of  Equation  (2.1)  for  breakdown  studies.  To  demonstrate 
this  fact,  let  us  consider  the  steady-state  equation  since  it  will 
be  shown  in  Section  2.2  that  the  determination  of  breakdown  condi¬ 
tions  involves  the  steady-state  solutions  to  Equation  (2.1).  Let 
Q  be  the  ionization  rate  due  to  collisions  among  the  neutral 
particles.  Then  the  complete  electron  conservation  equation  in 
steady  state  becomes 

-  •  ?nx  +  V2  (Dn^ )  +  vnA  +  Q  -  0  (2.2) 

where  we  have  used  the  abbreviation  v  ■  (Vj  -  vA  -  vR)  and  n^ 
is  the  electron  density. 

Consider  now  the  electron  conservation  equation  in  the  absence 
of  the  applied  electric  field  -  for  example,  that  which  could  apply 
far  upstream  of  the  electric  field,  namely, 

-  (V  •  V)n0  +  (D0n0)  ♦  (vt  -  vA  -  vE)o  nQ  +  Q  ■  0  (2.3) 

The  terms  Q  in  these  two  equations  are  assumed  equal  because  each 
depends  on  the  energy  of  the  heavy  particles,  an  energy  which  is 
assumed  unaffected  by  the  presence  of  the  field.  The  applied  field 
heats  the  electrons  largely,  not  the  ions,  and  the  electrons  cannot 
readily  transfer  their  energy  to  the  heavy  particles  (Refs.  1,  12, 
and  13).  If  Equation  (2.3)  is  subtracted  from  Equation  (2.2),  one 
obtains 

D  ^  n  +  vn  + 

o)  oj 

where  n  =  n  -  . 

1  o 


-  (V  •  V)n  +  V2 


Dn  +  (r 
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Within  the  applied  field,  the  solution  we  seek  (the  unique 
solution  between  breakdown  and  squelchinq  discussed  in  Section  2.2) 
will  generally  have1  n  >>  nQ.  On  the  other  hand,  just  outside 
the  applied  field  we  have  v  »  v  -  v  ,  D  »  (D  -  Dq) ,  with 
n  _>  nQ.  Therefore,  in  all  regions  of  interest,  the  last  equation 
may  be  written  to  sufficient  accuracy. 

-  (V  •  V) n  +  V2 (Dn)  +  vn  -  0  (2.4) 


This  is  the  same  as  Equation  (2.1)  in  the  steady  state.  In  addi¬ 
tion,  we  see  that  since  the  boundary  conditions  for  Equation  (2.2) 
are  n  *•  nQ  at  x  ■  +_  »  or  z  ■  +  <*>;  the  boundary  conditions 
for  Equation  (2.4)  are  n  ■  0  at  x  ■  ;*  <»  or  z  "  +_  «>.  As  is 
common  in  breakdown  studies,  the  surface  is  assumed  to  absorb  (or 
annihilate  by  acting  as  a  third  body  in  recombination)  all  elec¬ 
trons  striking  it.  Thus,  we  take  n^  “0,  nQ  "  °>  and  n  “  0 
at  y  ■  0  as  well  as  at  y  ■  <». 

2.1.2  The  linear  electron  loss  term 

The  electron  attachment  is  linear  in  n  if  either  dissocia¬ 
tive  attachment  predominates  over  three-body  processes  or  if  the 
degree  of  ionization  is  so  low  that  heavy  particles  predominate 
over  electrons  as  third  bodies  during  attachments .  This  will  be 
true  when  x  "  VN  <  ~  10-4 ,  where  N  is  the  heavy  particle 
density.  The  recombination,  however,  will  be  proportional  to  n2 
and  even  n3  for  very  large  x* 

Experimentally  measured  values  of  attachment  coefficients 
in  cold  air  indicate2  (Ref.  2b,  p.  150) 
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_!  dn\ 


n  dt  / 


attach 


2xl04  p 


sec 


-l 


1See  Reference  8  for  a  typical  solution. 

2Here  Eg  is  the  "effective"  field  strength  discussed  in  Section  3.1.3 
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where  E_  is  in  volts  per  centimeter  and  p  is  in  mm  Hg,  whereas 
experimental  values  of  recombination  coefficients  in  cold  air  give 
approximately  (Ref.  2b,  p.  159) 


(k  drA  _  (l 00 

^d4ecomb  U  2‘3J 


xlO-6  n  sec 


with  n  in  cm-3.  Thus,  when  p  is  of  the  order  of  1  mm  Hg  and 
E  is  greater  than  10  volts  per  centimeter,  the  ratio  gives 


Attachments  .  „  Ee  <v°HVc») 

Recombinations  '  p(mm)n(cm-s) 

In  the  other  extreme,  for  E^p  small  with  p  large,  the  ratio 
gives 

Attachments  =  lO10  Pt™1*) 

Recombinations  ”  ,  -3. 

n  (cm  ) 


Thus,  at  least  for  cold  air,  recombinations  are  negligible  com¬ 
pared  with  attachments  except  when  the  degree  of  ionization  is 
very  high. 

In  addition,  for  the  reentry  problem,  we  are  largely  inter¬ 
ested  in  cases  where  diffusion  loss  predominates  so  that  the  effects 
of  both  vA  and  vR  on  electron  loss  will  be  small  compared  with 
diffusion  losses. 

2.1.3  The  form  of  the  diffusion  term 

The  diffusion  term  in  Equation  (2.1)  is  Vs (Dn)  rather  than 

(D?n)  as  commonly  encountered  in  gas  kinetics.  It  is  shown 

in  Sections  3.2.1  and  3.2.2  that,  when  (1)  the  collision  frequency 

between  electrons  and  neutral  particles  is  much  greater  than  the 

collision  frequency  between  electrons,  that  is,  v  »  v  and 

en  6e 

when  (2)  the  average  momentum  transfer  collision  frequency  changes 
slowly  with  position  in  the  region  of  interest,  the  form  of  the 
diffusion  term  is  indeed  V2 (Dn)  and  that  this  term  includes 
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thermal  diffusion.  (The  main  physical  reason  for  the  difference 
between  the  structure  of  the  electron  diffusion  term  and  the  usual 
gas  kinetic  diffusion  term  is  that  the  electron  mass  is  so  extremely 
small  that  it  does  not  contribute  appreciably  to  the  average  mass- 
flow  velocity.) 

Let  us  now  examine  whether  the  restrictions  just  mentioned 
are  fulfilled  in  the  reentry  application.  The  usual  estimate  for 
v„  (Ref.  14)  gives 

G© 


ee 


2?rn 


_2  _  3 

m  c 


In 


1  + 


2  _  4  f  kT 
me/  e 


.  4Trne2 


sec 


-l 


(2.5) 


where  all  of  the  symbols  are  defined  in  the  list  of  notations.  If 
we  take  a  cross  section  of  10“15  cm2  for  electron-neutral  colli¬ 
sions,  then 


v0n  =  N(10-15)  YikT^/irm 


and,  noting  that  c  «  V^e  » 


vee  _  10xoy 
*en 


(2.6) 


where  y  "  >VN.  Thus,  vee  «  ven  is  achieved  if  y/Te2  «  10-10. 
This  is  easily  fulfilled  for  temperatures  of  thousands  of  degrees, 
such  as  are  common  in  reentry. 

Prom  Equation  (3.31)  it  will  be  seen  that  the  second  require¬ 
ment  mentioned  above  is  approximately  equivalent  to 

~r\  (2.7a) 

c2  / 


W 


en 


en 


«  (  the  larger  of  ~ 


The  data  and  estimates  of  Reference  15  indicate  that  v _  is 

en 

fairly  independent  of  energy.  Thus,  the  condition  of  Equation  (2.7) 
is  largely  density  dependent  and  may  be  written 


I 


2-6 


^  «  ^.arger  of  or  ^  (2.7b) 

It  is  quite  possible  that  this  condition  is  not  always  satisfied 
in  a  hypersonic  boundary  layer  and  thus  the  term  omitted  from 
Equation  (3.31)  to  obtain  the  form  ^Dn  in  Equation  (2.1)  should 
be  included.  However,  the  data  and  estimates  available  for  D, 
ven>  and  c2  are  so  poor  under  conditions  obtained  in  a  reentry 
boundary  layer  that  such  a  correction  to  our  breakdown  calculations 
would  have  little  meaning.  Since  the  omitted  term  is  linear  in  n, 
there  would  be  no  problem  in  including  such  a  term  in  the  procedure 
developed  in  Section  5  for  the  computation  of  breakdown  field 
strength,  that  is,  if  the  needed  data  were  available. 

For  the  reentry  problem,  the  diffusion  will  quite  often  be 
fully  ambipolar.  For,  as  shown  in  Section  3.2.3,  if  n  is  greater 
than  3xl03  era"3,  the  ambipolar  diffusion  is  fully  established. 

Such  a  minute  electron  density  is  certainly  attained  in  many  reentry 
boundary  layers.  However,  to  check  our  computational  methods  with 
experiment,  it  was  necessary  in  Section  6  to  make  a  few  computations 
using  the  free-dif fusion  coefficient,  since  the  only  experiments 
involving  the  geometries  of  concern  here  were  for  cold  un-ionized 
air  under  no-flow  conditions  (Ref.  4). 

2 . 2  The  Breakdown  Concept 

2.2.1  Continuous -wave  breakdown  _ 

We  shall  first  discuss  breakdown  due  to  continuous-wave  opera¬ 
tion  of  the  transmitting  antenna.  In  Equation  (2.1),  both  macro¬ 
scopic  parameters  D,  the  diffusion  coefficient,  and 
v  5  Vj  -  vA  -  vR,  the  net  ionization  coefficient,  are  functions 
of  the  average  electron  energy.  We  will  assume  that  the  electron 
random-velocity  distribution  is  "relaxed"  or  "stationary"  so  that 
the  average  electron  energy  is  a  function  only  of  the  "effective" 
electric  field  strength  and  species  concentrations  at  any  position 
in  question.  The  discussion  of  Section  3.1.2  examines  the  conditions 
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under  which  this  assumption  is  valid.  It  is  found  that  a  large 
portion  of  the  reentry  regime  fulfills  this  condition.  We  also 
assume  that  the  mixture  of  heavy  particles  is  not  appreciably 
changed  by  the  applied  field  and  resulting  ionization.  Thus,  all 
coefficient  functions  in  Equation  (2.1)  are  functions  only  of 
position  and  the  effective  electric  field  strength. 

At  a  fixed  position,  Equation  (2.1)  behaves  in  a  manner 
similar  to  the  expression 


dn 

3t 


sn 


where  s  is  negative  for  weak  electric  fields  and  positive  for 
strong  electric  fields.  Thus,  there  are  only  three  types  of 
steady-state  solutions:  n  ■  0,  corresponding  to  s  negative 
and  no  breakdown;  n  ■  0,  corresponding  to  s  positive  and 
giving  breakdown;  and  a  unique  condition  at  which  s  ■  0,  or, 
from  Equation  (2.1) 

V2  (Dn)  -  V  •  Tn  -  -  vn  (2.8) 

This  equation  determines  the  electric  field  strength  above  which 
breakdown  will  occur.  Once  ionization  growth  sets  in  at  a  rate 
above  that  occurring  in  the  solution  to  this  equation,  the  gas 
cannot,  in  steady  state,  flow  through  the  applied  field  fast  enough 
to  limit  the  ionization  to  a  finite  value.  Ionization  would  remain 
finite,  but  the  gas  would  become  completely  ionized,  if  in  v  we 
had  accounted  for  the  depletion  of  the  number  of  ionizable  parti¬ 
cles.  The  ionization  becomes  infinite  everywhere  due  to  the  action 
of  the  diffusion  term  coupled  with  the  unrealistically  unlimited 
supply  of  ionizable  particles. 

Let  us  put  Equation  (2.4)  or  (2.8)  in  a  form  that  will  make 
use  of  tabulated  data  (for  cold  air)  more  convenient.  Take 
rf/  -  nD,  then,  since  in  the  problem  considered  here,  •  V  m  0, 
the  equation  may  be  rewritten  as 
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& 

c.> ■ - xs  (f Xv)  * <2-9> 

where  £  =  v/DE  8  is  the  high-frequency  ionization  coefficient 
for  the  mixture  in  question,  and 

k8  -  C'E  '8  (2.10) 

e 

Here  the  effective  electric  field  is  given  by 

Ea  -  E  'F(x,y,z)  (2.11) 

where  Ee'  is  the  effective  field  at  some  specified  point  in  the 
region  of  interest  and  £'  is  the  value  of  £  at  this  point. 
Determination  of  the  field  strength  above  which  breakdown  occurs 
then  involves  determining  the  lowest  eigenvalue  k2  of  Equation 
(2.9)  and  thus,  E  through  Equations  (2.10)  and  (2.11).  In  the 
case  of  a  high-frequency  alternating  electric  field,  Ee  is  the 
effective  d-c  field  strength  defined  in  Section  3.1.3. 

Let  us  use  the  abbreviations 

G-  (£/£')  (E@/Ee  * ) 2  and  B  -  ?/D  (2.12) 

so  that  Equation  (2.9)  is  simplified  to 

7 -  7  •  fify)  -  -  k8G^  (2.13) 

2.2.2  Multiple -pulse  breakdown 

As  in  the  c-w  case  we  are  again  concerned  with  the  condition 
that  ionization  does  not  exhibit  unlimited  growth  as  time  increases. 
Let  the  discussion  for  the  moment  be  restricted  such  that  D  is 
assumed  constant  in  time.  The  discussion  with  nonconstant  D  will 
follow. 


I 
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|  Suppose  we  assume,  as  in  Sections  2.2.1  and  2.2.2,  that  v 

is  a  function  only  of  the  effective  local  electric  field  and  the 
I*  local  mixture.  As  stated  earlier,  this  necessitates  the  assumption 

of  a  stationary  or  relaxed  random-velocity  distribution  function 
(RVDF) .  The  discussion  in  Section  3.1  will  show  that  for  pulses 
! .  this  assumption  places  certain  restrictions  on  the  pulse  width  and 

on  the  field  strength.  Let  us  expand  v  in  a  series 

v "  Z  ri(^)EeYi  (2.14) 

i 

and  let  the  effective  electric  d-c  field  (defined  in  Sections 
3.1.3  and  3.2.4)  be 


Ee  -  Ee'F(r)P(t)  (2.15) 

where  Ee'  is  the  value  of  Eg  at  some  point  7*  and  P(t) 

gives  the  repeated  pulse  shape.  For  a  square  pulse,  if  e(t)  is 

a  unit  step  function,  the  function  P(t)  is 

P(t)  -  £e(t)  -  s(t  -  T)  +  e(t  -  T)  -  e(t  -  T  -  t)  +  ... 

+  e(t  -  qT)  -  e(t  -  qT  -  T)  +  ...J  (2.16) 

where  q  is  an  integer,  t  is  the  pulse  width,  and  T  is  the 
pulse  spacing.  Note  that  if  P  is  raised  to  any  power,  then 

p7  -  p 

because  of  the  nature  of  step  functions . 

Equation  (2.1)  may  be  rewritten  as  (D  is  assumed  constant) 


+  V  •  Vn  -  DV2n  -  (k*)*G  (7,7' ) P(t) n 


(2.17) 
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whs  re 


and 


G  (7,7' ) 


r  t  v  * 

]jTr.(7)  e6'f(7)  1 

£>*(?•)  (Be')n 


(2.18a) 


(2.18b) 


Consider  now  the  Laplace  transform  7  of  any  function  f(t) 
that  is  either  periodic,  or  constant  and  finite,  or  diverges  more 

rv  4* 

slowly  than  e  for  large  time.  Then 


00 

?-J  e“St  f(t)dt  ,  Re (s)  >0  (2.19) 

o 


This  definition  may  be  used  even  for  f  diverging  as  fast  or 

rv  f 

faster  than  e  for  large  t,  so  long  as  we  do  not  attempt  to 
invert  the  transform.  Then,  of  course,  ?  ■  <».  For  such  a  func¬ 
tion,  the  following  theorem  may  be  proved: 


lim 
s  -*  0 


lim 
t  —  00 


r  t+T 

If  J  f(t)dt 

t 


if  f(t)  goes  to  a  "n 
constant,  is  periodic 
with  period  T,  or 
diverges  more  slowly 


than 
t 


-at 


for  large 


>(2.20) 


V- 


00 


if  f(t)  diverges  as 
fa|t  or  faster  than 


Let  us  call 


[•*«•)] 


(2.21) 
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Take  the  Laplace  transform  of  Equation  (2.17),  obtaining 
-  nt-Q  +  ^  •  Vn  -  D72n  -  k*2G  p(t)n 


sn 


Multiply  through  by  s  and  take  the  limit  as  s  -*  0.  Then  if  f 
does  not  diverge  exponentially,  using  Equations  (2.20)  and  (2.21), 
one  obtains 


V  •  ?  ^ny  -  D72  (n^  -  k*2G  n 


("*) 

<n> 


(2.22) 


If  n  is  periodic  at  large  t,  it  must  have  the  same  period  as 
P,  namely,  T,  such  that 


<">  ■  T  / 

qT 


(qfi)T 


n  (t)dt 


(2.23a) 


and 


qT+T 


<*”>  -  «1- -  ?  f  » 

qT 


(t)dt 


(2.23b) 


since  P(t)  vanishes  for  qT+T<t<(q+l)T.  At  fixed  T 
within  the  field,  while  the  field  is  on,  n  behaves  something  like 

s  (t-qT) 


n  ~  n  e 
o 


qT  ^  t  £  (qt  +  T) 


(2.24a) 


and  while  the  field  is  off,  it  behaves  like 

-s  (t-qT-x) 


n  ~  n 


(t-T)e 


,  (qT  +  T)  ^  t  £  (q  +  1)T  (2.24b) 


But  at  critical  conditions  between  breakdown  and  no  breakdown,  one 
must  have 


n(t-qT)  "  n  t-(q+i)T 


(2.24c) 
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as  depicted  schematically  in  Figure  2.2.  Then  from  these  last 
three  relations,  one  obtains 

st  -s  (T-t) 

~  °  _  1  _  i 


or 


s 


o 


(2 . 24d) 


Then  using  (2.24a)  and  (2.24b)  and  (2.23a)  and  (2.23b),  one  finds 


Then  with  these  relations  and  Equation  (2.24d),  one  finds 


Substitution  of  this  into  Equation  (2.22)  gives 

(V  •  V)  (n^  -  D  V2  (n)  -  k*2  ^  G  (n^  (2.25) 

The  eigenvalues  of  this  equation  give  the  breakdown  field 
strength  just  as  they  do  for  Equation  (2.8)  or  (2.13)  in  the  c-w 
case.  In  addition,  for  constant  D,  the  equations  are  identical 
except  for  the  factor  t/t  on  the  right-hand  side.  Notice  that 
k*2  defined  above  is  Dk2.  Thus  the  multiple-pulse  c&se  is  iden¬ 
tical  to  a  c-w  breakdown  problem  with  k*2G  ■  v  reduced  by  the 
factor  t/T. 
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For  the  case  of  nonconstant  D,  the  treatment  is  very  simi¬ 
lar.  The  transformed  conservation  equation,  multiplied  by  s, 
becomes  in  the  limit  as  s  0 


(V 


?>  <-)  - 


<&\ . 


Jt*sG  (n) 


(n) 


(2.26) 


where  ^Pn^/^n^  is  again  t/t.  The  new  effective  diffusion 
constant 


eff 


<»> 


(2.27) 


may  be  evaluated  in  a  manner  similar  to  ^Pn^/^n^.  For  the  relaxed 
RVDF  (Section  3.1),  D  is  a  function  only  of  the  average  electron 
energy,  and  the  average  electron  energy  ia  determined  solely  by  the 
effective  applied  field  strength.  Therefore,  for  a  square  pulse  at 
a  fixed  position  in  space,  D  is  discontinuous  in  time;  it  has  one 
value  when  the  field  is  on  and  another  when  the  field  is  off.  The 
size  of  the  discontinuity  depends  upon  position  because  of  the 
antenna  pattern.  Such  an  evaluation  is  then  straightforward  but 
cannot  be  carried  out  very  precisely  since  the  dependence  of  the 
average  electron  energy  on  E  is  not  very  well  known.  In  addition, 
there  i3  the  problem  that,  when  the  field  is  off,  the  electrons  do 
not  cool  very  readily  back  to  the  average  field-off  energy;  the 
cooling  time  is  approximately 


(H/2m)  (l/3ven) 


(2.28) 


where  M  is  the  heavy  particle  mass,  m  is  the  electron  mass, 

and  v  is  the  electron-neutral  collision  frequency.  For  inter- 
en 

vals  between  pulses  shorter  than  this  time,  the  value  of  D  with 
the  field  on  and  constant  in  time  would  probably  be  a  reasonable 
approximation.  Then  one  can  take 


<d„> 


=  D  (field  on) 


(2.29) 
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and  again  we  have  the  result  that  the  multiple-pulse  problem  is 
identical  to  a  c-w  eigenvalue  breakdown  problem  with  an  effective 
net  ionization  frequency  of  (t/t)v .  This  result  was  obtained  in 
Reference  3  with  what  appears  to  be  a  much  simpler  argument; 
however,  it  is  believed  that  the  present  treatment  brings  out  the 
validity  of  the  corrected  eigenvalue  approach  more  clearly. 

It  then  appears  that,  as  the  pulse  interval  goes  to  infinity 
(T  -*  oo)  ,  the  eigenvalue  and  the  electric  field  for  breakdown  goes 
to  infinity.  In  such  a  case,  however,  one  must  consider  that  the 
problem  has  gone  over  into  single-pulse  breakdown  and  thus  the 
breakdown  field  strength  is  limited  to  the  value  determined  there 
(see  the  next  section) .  There  would  have  been  no  need  for  special 
treatment  of  the  single-pulse  case,  that  is,  it  could  be  treated 
naturally  in  the  limit  as  T  becomes  large,  if  in  the  treatment 
of  the  random-velocity  distribution  function  (Section  3),  and  thus 
in  v ,  we  were  to  include  the  effect  of  plasma  resonance  in  the 
rate  of  heating.  However,  to  include  plasma  resonance  in  an  accu¬ 
rate  way  in  the  kinetic  equation  introduces  a  considerably  more 
complex  kinetic  equation  (see,  e.g.,  Appendix  A). 

2.2.3  Single-pulse  breakdown 

Since  it  is  difficult  to  include  plasma  resonance  in  the 
kinetic  equation  for  evaluation  of  the  rate  of  heating,  because 
of  the  nonlinearity  introduced,  the  concept  of  single-pulse  break¬ 
down  has  been  used  in  the  literature  (Refs.  3,  9,  and  11).  Here 
one  is  no  longer  concerned  with  an  eigenvalue  problem,  but  rather 
with  limitations  imposed  by  the  condition  that  resonant  plasma 
oscillations  do  not  occur.  If  such  oscillations  were  to  occur, 
resonant  absorption  of  the  applied  energy  by  the  electrons  would 
result  (Ref.  16,  p.  50),  and  thus  the  ionization  growth  rate 
would  be  greatly  increased. 

The  condition  that  such  resonance  does  not  occur  is  that  the 
applied  frequency  be  greater  than  the  plasma  frequency 


go 

2tt 
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2v 


(2.30a) 
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where 

-  8.97X103  npl/s  sec-1  (2.30b) 

Here  n  is  the  density  of  electrons  (in  cm-3)  giving  rise  to  the 

Sr 

plasma  frequency  o>p. 

The  breakdown  field  strength  is  evaluated  from  this  condition 
by  using  Equation  (2.1)  in  a  coordinate  system  with  zero  velocity, 
that  is,  moving  with  the  flow.  The  condition  that  no  parcel  of 
air  should  attain  the  electron  density  np  during  the  period  T 
of  the  applied  field  becomes,  by  integration  of  Equation  (2.1) 
without  the  velocity  term, 

*”  (%)  *  /  [VlEe»  +  2f^2ai]dt  (2'31) 

where  nQ  is  the  initial  electron  density.  To  evaluate  this  rela¬ 
tion,  one  takes  72(Dn)/(Dn)  ■  1/AZ,  where  A  is  the  diffusion 
length  which  may  be  evaluated  from  solution  of  the  steady-state 
problem  (Eq.  (2.8)),  or  approximated  by  the  length  over  which 
escaping  electrons  must  diffuse  to  get  out  of  the  applied  field. 

To  compensate  for  the  neglect  of  resonant  plasma  oscillations 
in  the  multiple-pulse  problem,  the  value  of  Eg  determined  from 
Equation  (2.31)  is  used  as  an  upper  limit  for  large  spacing  between 
pulses.  Typical  multiple-pulse  results  using  this  upper  limit, 
computed  in  Reference  3,  are  shown  in  Figure  2.3.  These  data  have 
been  obtained  from  breakdown  measurements  and  theory  in  a  closed 
cavity. 
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Figure  2.1.-  Coordinates  and  geometry. 


Critical  condition 


No  breakdown 


Figure  2.2.-  Multiple-pulse  breakdown  schematic. 
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3.  ESTIMATES  OP  MACROSCOPIC  PROPERTY  INPUTS  FROM  THEORY  AND 

EXPERIMENT 

In  this  section,  the  theoretical  relationships  between  micro¬ 
scopic  collisions  and  macroscopic  properties  is  reviewed  and  used 
to  help  in  estimating  macroscopic  properties.  It  is  also  used  to 
determine  restrictions  on  the  applicability  of  experimental  data 
obtained  with  various  laboratory  experiments  and  environments .  In 
particular,  it  is  pointed  out  that  the  usual  d-c  laboratory  data 
have  all  been  obtained  under  conditions  of  a  relaxed  electron 
random-velocity  distribution  function  (RVDF)  and  are  applicable  only 
under  conditions  of  a  relaxed  RVDF.  Also  pointed  out  is  the  impor¬ 
tance  of  using  data  for  the  precise  (or  nearly  precise)  mixture  in 
question. 

The  criteria  for  this  relaxed  RVDF  are  analyzed,  and  it  is 
demonstrated  that  many,  although  not  all,  of  the  reentry  situations 
involve  relaxed  distributions.  For  the  case  of  unrelaxed  distribu¬ 
tions,  it  is  pointed  out  that  the  breakdown  problem  is  still  com¬ 
pletely  tractable  with  the  use  of  reasonable  functions  for  the 
microscopic  ionization  frequency,  excitation  collision  frequency 
and  attachment  frequency,  as  well  as  the  microscopic  momentum  trans¬ 
fer  collision  frequency.  These  microscopic  quantities  are  functions 
of  the  individual  electron  velocity;  that  is,  they  are  not  yet  aver¬ 
aged  over  a  velocity  distribution.  It  is  demonstrated  that  the 
unrelaxed  breakdown  problem  may  be  handled  with  a  variational  prin¬ 
ciple  similar  to  that  used  in  this  report. 

It  is  believed  that  the  above  points  have  not  been  previously 
discussed  in  the  literature.  In  addition,  some  theoretical  rela¬ 
tions  have  been  misinterpreted  or  misused  in  portions  of  the 
literature.  For  clarity,  therefore,  results  have  been  rederived 
here.  In  particular,  one  misconception  is  that  0fehe  electron  RVDF 
is  either  "Maxwellian"  or  "Druyvestein"  or  somewhere  in  between 
these .  Another  misconception  is  that  in  the  relationship 
D  "  2/3  p.u,  between  the  diffusion  constant,  the  mobility,  and 
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the  average  particle  kinetic  energy,  the  factor  of  two-thirds  is 
exact  only  for  a  Maxwellian  distribution. 

3.1  Kinetic  Equations  and  Velocity  Distribution  and  Density 

Functions 

3.1.1  The  reduced  kinetic  equation  and  the  random-velocity 

function  (RVDF) 

For  prediction  of  macroscopic  properties  from  microscopic 
collision  phenomena,  one  must  know  the  velocity  distribution  of 
the  electrons.  The  form  of  the  velocity  distribution  must  be 
determined  from  the  kinetic  equation.  In  Reference  1,  12,  13,  14, 
17,  18,  or  19,  it  is  shown3  from  the  kinetic  equation  that,  under 
certain  conditions,  the  electron  density  function  in  velocity  and 
position  space  is  given  approximately  by 


F  (t,o),c,r) 
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*  Fo°(t,u),  |'c|  ,7) 
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en 
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en 
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(3.1) 


where  c  is  not  the  actual  electron  velocity,  but  rather  the 
random  velocity  defined  in  terms  of  the  actual  electron  velocity 
7  and  the  macroscopic  flow  velocity  7  by 

7  -  7  -  7  (3.2) 

and  the  electric  field  is  assumed  to  be  E  e3a>t.  The  function 

o 

Fq°  satisfies  the  reduced  kinetic  equation  (assuming  7  ■  Vic) 


+  7  •  VF  ° 
o 


en 


+  or 


kT 

me 


(3.3) 


3In  Appendix  A,  some  of  the  special  restrictions  and  considerations 
imposed  on  this  kinetic  development  by  the  flow  and  geometry  in  the 
present  application  are  considered,  as  well  as  a  review  of  the 
development. 
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The  terms  on  the  right  are  due  to  collisions;  the  first  is  due  to 

elastic  collisions  and  the  second  is  due  to  all  inelastic  collisions 

(and  reappearance  after  an  inelastic  collision) .  Here  v^,  vA, 

vD,  and  v  _  are  the  ionization  frequency,  etc.,  for  an  individual 
k  en 

electron  cf  velocity  c,  as  opposed  to  the  average  frequencies 
(with  the  bars)  used  earlier.  The  function  q(c)  is  the  appearance 
or  reappearance  rate  for  electrons  that  have  been  produced  or  lost 
energy  in  an  inelastic  collision.  The  density  function  Fq°  is 
normalized  such  that 

Fq°  dcx  dcy  dc2  -  n(7,o),t)  (3.4) 

The  most  restricting  of  the  conditions3  under  which  Equations 
(3.1)  and  (3.3)  are  valid  are  (i)  the  requirement  of  low  ionization 
density  such  that  electron-electron  collisions  are  neglected,  and  thus 


(3.5) 


and  (ii)  a  high-frequency  electric  field  such  that  relaxation  is 
negligible  within  a  cycle,  so  that 

?en'  uhere  (*«  '  %n  *  *.l)  <3-6> 

The  restrictions  imposed  by  Equation  (3.5),  when  evaluated  in 
the  same  manner  as  (2.6)  with  (2.5),  give  approximately 


X  «  10"14  Te2(°K) 


(3.7) 


where  T@  is  representative  of  the  average  electron  random  energy 
such  that  (l/2)mcs  =  3'kTe/2.  Since  Tg  will  usually  be  between 
about  2,000°  and  15,000°  K  (as  will  be  discussed  later),  this 


3See  footnote  3  on  page  3-2. 
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restricts  the  model  to  a  degree  of  ionization  significantly  below 
10”5  to  10”8.  Such  a  restriction  is  probably  fulfilled  in  reentry 
along  the  sides  of  the  body,  that  is,  if  the  electrons  are  not 
"frozen  in"  at  a  level  corresponding  to  stagnation  temperature. 

At  stagnation  conditions  the  reentry  degree  of  ionization  may  reach 
and  even  slightly  exceed  these  values.  If  this  happens,  but  the 
ionization  is  still  low  enough  to  satisfy  vee  «  ven  such  that 

10"10  T  2  (°K)  »  X  ~  10“14  T  *{°K)  (3.8) 


then,  one  obtains  to  sufficiently  good  approximation  (Refs.  17  and 
18)  an  equation  identical  to  Equation  (3.3)  except  that  the  mass 
of  the  heavy  particles  is  replaced  by 


M* 


(3.9a) 


and  the  temperature  of  the  heavy  particles  is  replaced  by 


T*  -  T  *  Te  5  (^)  “  (nSf)  (3-9b> 

where 

«(?)  -  erf  *  -  2L-.exp(-  £2)  (3.9c) 
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Here,  as  before,  3kTe  =  me2.  The  equation  is  actually  nonlinear, 
since  one  must  use  Fq°  to  evaluate  c2;  however,  one  can  approxi¬ 
mate  the  c2  expected,  solve,  and,  if  necessary,  iterate. 
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The  function  ft(£)  is  very  small  for  £  <  1  and  grows  to  a 
maximum  of  unity  when  £  is  very  large.  The  effect  of  the 
electron-electron  collisions  is,  therefore,  to  make  the  average- 
and  faster-than-average-velocity  electrons  act  as  if  they  see 
lighter  but  hotter  neutrals.  Since  an  electron  could  transfer 
energy  more  easily  to  lighter  and  with  more  difficulty  to  hotter 
particles  by  collision,  these  two  effects  tend  to  compensate  one 
another.  This  possibly  helps  to  explain  why  breakdown  predictions 
have  been  so  accurate  even  with  neglect  of  electron-electron  col¬ 
lisions  in  the  analysis;  however,  the  reduction  in  mass  effect 
surely  predominates . 

In  solving  these  equations  it  has  become  common  to  define  a 
random-velocity  distribution,  f,  given  by 

PQ°  -  fn (r , t)  (3.10) 


If 

average 


T)  is  any  function  of  the  magnitude 
rj  is  given  by 


j  "o’  |  =  c ,  then  the 
(c)Fo°(c)4irc2dc 

(3.11) 


Multiply  Equation  (3.3)  by  47rc2dc  and  integrate  from  zero  to  <». 
Then  with  use  of  Equations  (3.4)  and  (3.11)  integration  by  parts 
gives 


+  •  ?n  -  7*  [(£)  »]♦▼•[«(£*£) 

■(vi-vA-vR)n  (3.12) 
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since  q  ■  2v ^  +  vx»  due  to  the  fact  that  for  each  fast  electron 
lost  in  ionization  there  reappear  two  slow  electrons,  whereas  for 
each  fast  electron  lost  by  excitation  there  reappears  one  slow 
electron.  The  terms  involving  d/dc  vanished  at  the  limits  of 
integration  because  at  c  ■  »,  p  °(»)  ■  0.  Thus,  we  have  derived 
the  electron  conservation  equation  except  for  one  extra  term  which 
will  be  discussed  in  Section  3.2.2. 

Now,  substituting  Equation  (3.10)  into  Equation  (3.3),  one 
obtains  with  use  of  Equation  (3.12), 


If  +  V-Vfiv.  &£).(vi+Va+Vr+  VJ[+VI-»A-;R-q)f 

("<=  if + 1?  f)]  (3'13a) 


2  _d_ 
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and 


-$  =  c2  . 

0  3ven  fn 


Cc.V.31m).b1  , 

n  v  3  ven  / 


(3.13c) 


Equation  (3.13a)  must  be  solved  in  order  to  determine  the  random- 
velocity  distribution  function. 

3.1.2  Criteria  for  a  relaxed  or  stationary  RVDF  and  difficulties 
introduced  if  not  applicable 

Virtually  all  work  on  breakdown  in  gases  which  has  utilized 
a  velocity-distribution  function  for  electrons  to  compute  macro¬ 
scopic  parameters  has  been  concerned  with  the  breakdown  condition 
under  nontransient,  nonflow  conditions.  Previous  work  has  there¬ 
fore  not  considered  velocity-distribution  transients  either  in 
time  or  position  space  (such  as  would  arise  from  flow  into  a  field) . 
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The  RVDF  was  determined  solely  by  the  balance  between  energy  input 
from  the  electromagnetic  field  and  energy  losses  due  to  elastic 
and  inelastic  collisions  and  diffusion  (the  more  energetic  parti¬ 
cles  diffusing  faster) .  The  equation  satisfied  by  such  an  RVDF 
is  Equation  (3.13a)  without  the  b/bt  and  the  "v  •  V  terms, 
namely , 


i  V  •  (<*>nf)  + 
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(3.13d) 


Before  attempting  to  estimate  under  what  conditions  such  a 
simplification  is  possible,  let  us  consider  the  complications  intro¬ 
duced  when  this  simplification  is  not  possible.  With  a  nonstation¬ 
ary  RVDF,  one  might  as  well  go  back  to  Equation  (3.3)  for  breakdown 
studies.  In  the  c-w  case,  the  5FQ0/dt  term  may  be  omitted  for 
determination  of  the  breakdown  condition  for  the  same  reasons  as 
discussed  in  Section  2.2  (dn/dt“  0  condition).  In  the  multiple- 
pulse  case,  one  may  use  the  time  transform  technique,  as  in 
Section  2.2.2,  to  obtain  an  equation  for  ^F0°^  without  the  d/dt 
term.  For  the  nonstationary  RVDF  the  breakdown  problem  again 
reduces  to  an  eigenvalue  problem,  but  now  it  is  further  complicated 
by  the  presence  of  an  additional  dimension  or  independent  variable, 
the  quantity  c,  and  differential  operators  involving  this  addi¬ 
tional  variable  that  are  not  self-adjoint.  This  problem,  however, 
is  still  quite  tractable  and  may  be  handled  by  the  variational 
technique  described  in  Appendix  B.  The  main  difference  between 
this  and  the  treatment  of  Section  5.2,  which  will  be  used  in  the 
present  study,  is  that  the  computations  involved  in  such  a  problem 
would  require  the  machine  storage  of  at  least  2sn+1  pieces  of 
information  if  there  were  sn  pieces  stored  for  the  relaxed 
n-dimensional  self-adjoint  RVDF  problem,  and  if  about  the  same 
accuracy  were  desired. 
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Although  this  extra  difficulty  is  far  from  prohibitive,  it.  is 
entirely  unwarranted  under  the  present  state  of  our  information 
(especially  if  the  criteria  for  a  stationary  RVDF  are  met) .  This 
is  because  of  the  very  poor  information  presently  available  for 
VI (c) ,  vA(c) ,  v  (c),  and  even  ven(c).  Also,  data  for  the 
relaxed  problem  are  more  readily  available  since  any  measurements 
of  macroscopic  properties  under  stationary  conditions  may  be  used 
again  for  macroscopic  properties  under  similar  stationary  condi¬ 
tions  (but  not  under  transient  RVDF  conditions) .  In  addition, 
fairly  good  values  of  the  macroscopic  parameters  D  and  v  can 
sometimes  be  estimated  without  precise  knowledge  of  the  RVDF. 

Let  us  now  attempt  to  estimate  under  what  conditions  this 
stationary  assumption  is  possible.  This  criterion  must  require 
that  the  balance  between  energy  input  from  the  field  and  energy 
lost  must  be  established  during  flow  through  the  field  in  a  neg¬ 
ligible  distance  compared  to  externally  applied  field  width  (in 
the  flow  direction)  and  in  a  time  short  compared  to  the  pulse  width 
(when  there  is  one) . 

If  we  assume  that  the  initial  distribution,  upon  entering  the 
field,  is  Maxwellian  at  temperature,  T,  and  if  this  is  substituted 
into  Equation  (3.13a),  there  results 


Here  we  have  neglected  inelastic  collisions  for  we  are  interested 
in  evaluating  this  expression  at  energies  less  than  the  threshold 
energy  of  the  lowest  excited  state.  Also,  diffusion  loss  was  neg¬ 
lected,  for  there  will  be  little  such  loss,  initially. 

This  may  be  also  rewritten  as 

C  dz  )initial  "  nkTV  °RErms2  ("  1  +  T  kT~ )  =  H  (3,14) 
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defining  H  and  oR  (later  shown  to  be  the  electrical  conduc¬ 
tivity)  eo  that,  initially, 


f  “  fz-o(c)e 


z/H  (c) 


where  H  is  the  initial  relaxation  length  and  we  have  taken  the 
field  to  start  at  z  ■  0.  For  the  pulsed  nonflow  case,  we  would 
have  obtained 


d  [n  f  _  V 
dt  H 


or 


We  see  that,  initially,  from  the  behavior  of  H,  the  dis¬ 
tribution  function  is  "pivoting"  about  the  average  energy;  it  is 
decreasing  below  and  increasing  above  the  average  with  time  or 
length  of  flow  into  the  field.  In  the  region  of  (1/2) me2  ■  ux, 
the  lowest  inelastic  energy  level,  the  requirement  that  |h|  «  Lg, 
where  L_,  is  the  length  of  the  field-containing  region,  gives 

Ci 


nu. 


o_E 
R  rms 


(3.15a) 


whereas  at  the  low  energies  this  becomes 

L. 


nkT 


0_E 
R  rms 


E 

«  — 


(3.15b) 


Conditions  common  in  the  reentry  problem  are  ■  2  cm, 
ux  ~  1  e.v.,  T  ~  2000°  K,  p  ■  1  mm  Hg,  V  =  105  cxo/aec,  and 
we  will  take 

v  ■  6xlO®(p/T)c  (3.16) 

en 

where  p  is  in  mm  Hg  and  T  is  in  °K.  In  Equation  (3.16)  one 
must  take  c  ■  for  Equation  (3.15a)  and  c  ~  Y  kT/m  for 
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Equation  (3.15b).  Substitution  of  these  numbers  gives  from 
Equation  (3.15a) 

Erms  >>  0,7  volt/cm  (3.17a) 

and  from  Equation  (3.15b) 

E».m«  »~  1  . 5  volts/cm  (3.17b) 

riuo 


Since  these  minimum  field  strengths  are  only  about  1  to 
2  percent  of  the  field  strengths  required  for  breakdown,  one  might 
infer  from  this  that  the  relaxed  distribution  is  fairly  well 
achieved  for  breakdown  problems  in  reentry.  One  should  note,  how¬ 
ever,  that  if  we  set  Lg/V  ■  t  and  take  t  to  be  1  microsecond, 
the  analogous  results  become 


and 


E  »  7  volts/cm 
rms 

Erms  »  volts/cm 


(3.18a) 


(3.18b) 


for  the  high-  and  low-energy  regions .  Thus ,  for  microsecond  pulses 
under  these  conditions,  the  stationary  RVDF  may  not  be  quite 
achieved. 

Although  this  estimate  only  includes  initial  tendencies,  it 
should  give  a  reasonable  estimate  of  the  relaxation  requirements. 

3.1.3  The  stationary  distribution  as  a  condition  for  use  of 
d-c  data 

The  use  of  the  macroscopic  Equations  (2.1)  or,  equivalently, 
(3.12)  does  not  in  itself  assume  a  relaxed  RVDF.  The  only  place 
the  RVDF  comes  into  such  computations  with  the  macroscopic  equa¬ 
tions  is  through  D  and  v ,  the  diffusion  and  net  ionization 
coefficients.  The  assumption  of  a  relaxed  RVDF  comes  in  indirectly 
then  when  one  uses  data  from  experiments  conducted  under  relaxed 
conditions . 
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Despite  restriction  (3.6),  it  is  found  in  Appendix  A  that  at 

zero  frequency,  to  -  o,  one  obtains  Equations  (3.1)  and  (3.3) 

again  with  to  -  o  and  E _ _  replaced  with  E,  .  Then,  it  is 

3  rms  a-c 

immediately  obvious  that  one  obtains  the  same  distribution  function 
and  macroscopic  coefficients  with  a  high-frequency  a-c  field  as 
with  the  equivalently  effective  d-c  field  given  by4 


(3.19a) 


Therefore,  for  cold  air,  one  can  use  d-c  measured  data  for  high- 
frequency  breakdown  prediction.  However,  care  must  be  taken  because 
all  existing  d-c  data  were  obtained  under  relaxed  RVDF  conditions 
and  should  be  used  only  when  relaxed  RVDF  conditions  obtain. 

The  existing  measurements  of  pulsed  breakdown  (e.g.,  Ref.  3) 
are  also  conducted  under  relaxed  conditions.  This  may  be  seen  by 
noting 

Vrr,.'  5  V  (Sr-)  <3'19b» 

\  en ' 


Then,  the  conditions  (3.15a)  and  (3.15b)  become,  respectively, 


and 


7  (6x1°8»  b 


(6X108)  ±  V  2mu3 


(3.20a) 


(3.20b) 


4A  small  correction  to  this  equivalence  will  be  discussed  in 
Section  3.2.5. 


3-12 


Where  Equation  (3.16)  was  used  and  c  -  ~^2\i/n.  For  example,  the 
data  of  Reference  3  were  for 

Ee 

30  <  <  50  volts/cm-mm  Hg 

1  <  p  <  1000  mm  Hg 


X  ■  10,  10.7  cm 


T  =  Lg/V  =  1  M-  sec 


Then,  Equation  (3.20b)  gives  with  u2  =  (3/2 )kT  and  E^p  =  30 


=  1000  » 

pVt" 


Equation  (3.20a)  is  even  easier  to  satisfy. 

3.1.4  Determination  of  the  relaxed  distribution  function 

Equation  (3.13d)  has  been  solved  exactly  for  the  case  of  con¬ 
stant  v@n  in  terms  of  hypergeometric  functions.5  The  method  is, 
in  general,  very  complicated  and  tedious.  However,  the  main  fea¬ 
tures  of  the  distribution  function  may  be  obtained  by  the  approxi¬ 
mate  method  of  Allis  (Ref.  13),  Brown  (Ref.  1),  and  Smit  (Ref.  21), 
which  is  much  easier  to  use.  We  describe  this  below: 

The  usual  procedure  is  to  divide  the  problem  into  two  dif¬ 
ferent  regions.  The  first,  an  inelastic  region,  corresponds  to 
electron  energies  above  the  lowest  inelastic  threshold  u^,  that 
is,  all  u,  such  that 


u  >  u 

—  x 


Reference  20.  This  work  is  reviewed  in  Reference  1,  pp.  546-549. 
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Electron  energies  below  this  level  are  considered  in  the  elastic 
regions.  The  slope  of  the  logarithm  of  the  distribution  function 
in  both  these  regions  is  required  to  be  continuous.  We  consider 
first,  very  briefly,  the  inelastic  region. 

When  inelastic  collisions  are  present,  they  overpower  all 
other  collision  processes  because  of  the  very  large  energy  loss 
per  collision.  Also,  vA  and  vR  and,  of  course,  q  are  small 
at  large  energies.  The  term  involving  n/M,  which  is  due  to 
heavy  particle  recoil ,  and  the  diffusion  term  are  also  neglected 
in  this  region.  This  reduces  (3.13d)  to  an  equation  that  is  solved 
by  use  of  the  substitution  \ix/sf  ■  e“8.  The  reason  for  this 
change  of  variable  is  to  facilitate  the  requirement  of  good  accu¬ 
racy  at  ux,  the  lowest  inelastic  level,  and  u^,  the  ionization 
energy,  while  simultaneously  requiring  f (®)  -  0.  The  resulting 
first-order  equation  for  ds/du  is  solved  with  a  power  series  in 
1/u.  One  coefficient  in  this  power  series  is  reserved  for  matching 
at  ux  to  the  elastic  range  solution. 

In  the  elastic  range  Vj  and  vx  are  zero;  however,  qf 
must  be  considered.  As  noted  in  Section  3.1.1,  q  ■  2Vj  +  vx. 

Also,  q,  the  electron  reappearance  term,  is  appreciable  only  at 
low  velocities,  as  are  vA  and  v@n.  Then,  multiplying  Equation 
(3.13d)  by  47rc2dc  and  integrating  to  a  value  of  u  ■  (1/2) me2 
corresponding  to  moderate  velocities,  one  obtains 

(7I  +  7X  -  *)  -  -  X  ven°3  (uc  Is  +  TT  £)  (3-21> 

where  4>  is  the  diffusion  term  given  by 

c 

$>  s  J  n  V  ‘  ^nf)4lrcZdc  (3.22) 

o 

This  term  is  usually  small  in  the  elastic  region.  It  is  much  larger 
in  the  inelastic  region  but  is  usually  neglected  there  because  of 
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the  overpowering  effect  of  the  inelastic  collisions.  If  considered, 
it  is  generally  replaced  with  some  simplifying  approximation. 

The  solution  to  Equation  (3.21)  is 


where 


(3.23a) 


(3.23b) 


This  equation  is  valid  only  if  c,  and  thus  w,  is  large 
enough  to  be  above  the  values  at  which  q(c) ,  vA(c) ,  and  vR(c) 
contribute  appreciably.  Here,  wQ  is  the  value  of  w  when 
u  ■  uq,  and  uQ  is  the  average  value  of  u  above  ux  at  which 
inelastic  collisions  occur.  It  is  generally  slightly  above  ux> 

In  general,  uq  must  be  determined  by  matching  slopes  at  ux  to 
the  solution  in  the  inelastic  region. 

3.1.5  The  form  of  the  distribution  function  and  the  special  cases 
of  Maxwellian,  Druyvestein,  and  Margenau  distributions 

By  using  the  value  of  uc  from  Equation  (3.13b)  in  Equations 

(3.23a)  and  (3.23b),  the  form  of  the  distribution  function  and 

various  special  cases  may  be  predicted. 

When  the  first  inelastic  energy  level  is  so  high  and  w  (u  ) 

o  o 

so  large  (compared  to  average  w)  that  the  integral  in  Equation 
(3.23a)  is  dominated  by  the  upper  limit,  the  integral  is  approxi¬ 
mately  constant,  and  then  f  is 


f  a  (const) e“w 


(const) exp 


(3.24) 


I 

I 

I 

1 

1 

I 

I 

1 


3-15 


This  "Margenau  distribution"  was  first  derived  in  Reference  22. 
When 


or 


Me2 

3m2  v  2 
en 


«  kT 


e2E  2 


e  m 
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en 


a„E 
R  rms 


«  v 


en 


/2m  \_3 
V  M  J2 


kT 


}  (3.25) 


J 


this  reduces  to  the  ordinary  Maxwellian  (temperature)  distribution. 
The  condition  (3.25)  is  equivalent  to  the  statement  that  the  rate 
of  heat  input  is  negligible  compared  to  the  average  energy  exchanged 
between  electrons  and  heavy  particles. 

When  the  inequality  of  Equation  (3.25)  is  not  satisfied  and 
both  sides  are  about  equal,  one  obtains  again  a  Maxwellian  distri¬ 
bution  at  the  new  electron  temperature 


y  kT  +  o_E  J 
2  R  rms 


(3.26) 


if  ven  is  approximately  constant.  Essentially,  the  average 
energy  has  been  increased  by  the  energy  necessary  for  the  electrons 
to  lose  an  amount  during  collisions  with  heavy  particles  equal  to 
what  they  gain  between  collisions  from  the  electric  field. 

If  one  assumes  that  the  mean-free  path  (rather  than  the  col¬ 
lision  frequency)  is  approximately  independent  of  velocity,  then 
v  =  c/l,  and  if  (3.25)  is  satisfied,  one  obtains  the  Druyvestein 

CIl 

distribution  (Ref.  23) . 

When  uQ  is  not  very  much  larger  than  the  average  u,  one 
obtains  a  definite  distribution  by  requiring  continuous  slope  at 
ux  and  matching  (3.23)  to  the  inelastic  region.  The  form  of  the 
distribution  function  in  such  a  situation  will,  for  a  strong 
electric  field,  be  rather  flat  and  fall  off  suddenly  at  u  ■  ux. 


r. 


3-16 


The  slope  at  u  ■  ux  determines  the  number  of  particles  diffusing 
into  the  inelastic  region  and  thus  undergoes  inelastic  collisions . 
This  behavior  is  depicted  schematically  in  Figure  3.1.  A  careful 
treatment  of  the  computation  of  f  for  air  using  some  recent 
(though  not  yet  fully  accepted)  data  for  the  rates  of  vibrational 
excitation  by  electron  impact  may  be  found  in  Reference  24.  Care¬ 
ful  computations  with  noble  gases,  and  Heg  gas  (He  with  about 
1-percent  Hg)  are  in  Reference  20. 

3.2  Macroscopic  Properties  and  Macroscopic  Equations  from  Theory 

Although  the  very  poor  state  of  our  knowledge  of  ven^c)»  the 
momentum  transfer  collision  frequency,  vI(c),  the  net  ionization 
frequency,  etc.,  prevents  a  truly  reliable  prediction  of  macroscopic 
properties,  useful  formulas  for  extrapolating  and  interpreting 
experimental  data  may  be  obtained  from  theory. 

3.2.1  The  macroscopic  equations 

In  Section  3.1.1  we  derived  the  electron  conservation  equation 
from  the  kinetic  equation.  This,  in  addition,  gave  the  macroscopic 
parameters  for  this  equation.  By  taking  velocity  moments  of  the 
kinetic  equation  (3.3)  ,  all  the  other  macroscopic  equations  could 
be  likewise  derived. 

An  alternate  and  completely  equivalent  procedure  would  be  to 
write  down  the  macroscopic  equations  from  continuum  mechanics  and 
evaluate  the  necessary  macroscopic  parameters  by  use  of  Equation 
(3.1).  It  turns  out  that  the  equations  dealing  with  electrons  are 
not  always  independent  of  the  equations  governing  the  heavier  ions 
and,  therefore,  an  equation  for  heavy  ion  conservation  will  be 
included.  The  equations  of  interest  are  those  which  have  been 
already  averaged  over  the  very  short  time  2ir/co  of  the  high- 
frequency  oscillation.  The  time  dependence  remaining  in  the  equa¬ 
tions  is  for  the  much  slower  changes  such  as  growth  of  ionization. 

From  continuum  mechanics ,  the  particle  conservation  equations 


are 
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<“>)  ■  <7I  -  7a  -  ;r)  n 

Ir  -  7  ■(«&) '  T-i* I  -  v‘a  -  V_R> " 

'  '  av  ' 

where  eZa„  is  the  average  charge  per  ion.6  For  breakdown  studies, 

QV 

Z_„  will  be  very  close  to  unity.  The  average  ionization  vT  is 
a  function  of  the  precise  mixture  at  every  point  and  should  go  to 
zero  as  the  number  of  ionizable  particles  is  depleted.  This  latter 
property  is  not  important  in  the  computation  of  breakdown  condi¬ 
tions.  The  bar  denotes  averages  over  the  RVDF  and  the  ^notation 
denotes  time  averages  over  At  ■  2ir/u>. 

Let  us  write  for  the  diffusion  relation  for  electrons 

/T7\  /  _*  n) 

^c^v^-V- - ^ - enelD  (3.28a) 

and  for  the  heavy  ions 

-  DVN 

<3-28b> 


(3.27a) 

(3.27b) 


It  +  7  *  (n 


Here  we  have  assumed  special  forms  for  the  diffusion  which  will  be 
justified  in  the  next  section.  The  field  ED  is  due  to  the  charge 
separation  formed  by  the  different  diffusion  rates  of  the  ions  and 
the  electrons.  It  is  given  by  Poisson's  relation 

?  •(=<&)  -  -  ne  +  sI2ave  (3'29) 

where  eQ  is  the  free  space  dielectric  constant. 


6ln  Equation  (3.27b),  the  simplifying  assumption  is  made  of  an 

average  ion  of  charge  Z  e  and  that  ions  of  the  opposite  charge 
are  not  counted  (attachmlnt) .  This  serves  the  purpose  here  since 
this  equation  is  really  only  used  for  the  discussion  in  Section 
3.2.3  of  ambipolar  diffusion.  More  generally,  one  has  a  separate 
equation  for  each  ion,  and  all  the  individual  processes  must  be 
considered  on  the  right-hand  side  of  each  equation. 


•  -’•V-  .  . 
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An  energy  equation  is  not  needed  in  this  model  since  the 
assumption  of  a  relaxed  RVDF  determined  the  energy  as  a  function 
of  local  conditions.  The  energy  is  determined  by  the  local 
balance  between  electromagnetic  heating  and  collision  losses. 

3.2.2  Macroscopic  properties 

Equation  (3.11)  defined  the  average  of  a  microscopic  quantity 
when  it  is  a  function  only  of  the  magnitude  of  "c  and  not  its 
direction.  For  a  more  general  function,  one  that  is  dependent  on 
the  direction  of  c  as  well  as  its  magnitude,  one  must  use  F@ 
as  given  in  Equation  (3.1)  to  compute  averages.  This  requires 


n  tj  (c) 


I 


(3.30) 


It  is  easily  seen  that  Equation  (3.11)  is  a  special  case  of  (3.30). 

The  subsequent  average  over  (2t r/co)  gives  n(r}(c)^.  Because  of 
the  form  of  (3.1)  it  is  always  possible  to  reduce  averages  of  the 
type  (3.30)  to  averages  over  the  RVDF  as  expressed  by  Equation  (3.11). 
Let  us  proceed  in  this  way  to  derive  the  necessary  macroscopic  rela¬ 
tions.  In  deriving  these  relations  the  following  identity  is  used 


c  ■  A(c)dcx 


dc 


z 


A(c)4TTC2dC 


where  A  is  any  vector  that  may  be  a  function  of  the  magnitude 
of  "c.  This  result  follows  by  symmetry.  Some  terms  also  require 
integration  by  parts . 

In  order  to  obtain  an  expression  for  the  right-hand  side  of 
Equation  (3.28a),  another  term  must  be  added  to  Equation  (3.1)  to 
account  for  the  d-c  field  E^ .  This  term  is 

(-  e=  '  Vmvenc)  (dFo°/ac) 
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as  is  evident  from  the  derivation  from  the  treatment  of  d-c  fields 
in  Appendix  A.  (Alternately  one  could  use  the  existing  last  term 
in  Eq.  (3.1),  use  "e  ■  VT"?,,.  cos  cut  in  all  the  derivations, 
and  afterward  replace  Erin8/ri  +  (^2/ven2 ^  I  with  E^,  as 
required  by  Eq.  (3.19).)  With  this  extra  term  in  Equation  (3.1) 
we  find 


The  alternating  current  term  dropped  out  in  the  average  over  2tt/o>. 
The  first  term  of  this  expression  vanishes.  The  other  two  may  be 
rewritten,  by  use  of  the  identity  above,  as 


4Trc2dc 


Using  Fq°  ■  nf,  integrating  the  second  term  by  parts,  and 
manipulation  gives 


+ 


ne 

m 


f47TC2dC 


or,  by  use  of  Equation  (3.11), 
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Generally  the  terms  involving  derivatives  of  In  v  are  quite 
small  compared  to  the  others,  so  that  Equation  (3.31)  is  written 


n 


Comparison  with  Equation  (3.28a)  gives 


D 


e 


G£) and 


(3.32) 


(3.33) 


Analogous  expressions  for  the  right-hand  Bide  of  Equation 
(3.28b)  are  obtained  in  kinetic  theory  textbooks  (Ref.  25).  They 
give 


-  V 


kT 


MIvIn 
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Ze  -=r 
kT  ED 


(3.34) 


The  thermal  diffusion  ratio  k^  is  generally  small  because  the 
ion  mass  is  about  equal  to  the  average  particle  mass  -  this,  in 
fact,  is  the  essential  reason  for  the  different  form  of  the  dif¬ 
fusion  terms  in  Equations  (3.28a)  and  (3.28b).  Comparison  with 
Equation  (3.28b)  gives 


.  *T 

Vm 
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MIvIn 


and 


(3.35) 
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The  average 
course,  given  by 


electron  induced  ionization  frequency  is , 
Equation  (3.30)  as 


of 


a> 

2tt 


z-rr/ui  r  » 

/  <  I  /  [_aIi^c>Ni j  cf4lro‘dc 

i  L  1  CX1 


dt 


(3.36) 


where  (l/2)mcT.2  is  equal  to  the  ionization  potential  of  the  i ^ 
species,  N.  is  the  number  of  neutrals  of  the  i  species,  and 

1  i.V 

0^(0)  is  the  cross  section  for  ionization  of  the  i^  species. 

The  instantaneous  conductivity  of  the  gas  may  be  found  by 
computing  the_current  density,  using  Equations  (3.1)  and 

(3.30)  to  obtain  c,  which  gives 


oE  =  j  s  -  ne  c  =  Re  |  (oR  +  jo 


(dR  +  j°l)  . 


.  ne' 


Re 


(ven  +  1®)  + 


dv 


en 


[_(ven  +  J“) 

The  second  term  is  generally  negligible  so  that 


2  dc 


.j"*? 


dR  + 


^  ”  (v.n  4  j") 

The  instantaneous  heating  rate  is  7.  if  or 

•  if  ■  ne  "c  •  if  -  Eq2  ^or  cos2u)t  -  y  Oj  sin  2cot^ 


so  that 


(j  ■  - 


E  2on 
rms  R 


(3.37) 


(3.38) 


(3.39) 


Since  v_„  varies  slowly  with  electron  energy  (Ref.  15) , 
©n 

it  is  common  to  use  in  these  expressions 


3-22 


This  is  sometimes  called  the  constant  collision  frequency  approxi¬ 
mation.  Notice  that  with  this  approximation  Equation  (3.33) 
becomes 


De  *  cV3ven  ’  T  (1/"‘%n)  (?  “=*)  '  T  (t  "=*)  <3-40b> 


3.2.3  Ambipolar  diffusion  phenomena 

If  in  Equations  (3.27a),  (3.27b),  and  (3.29),  with  Equations 
(3.28a)  and  (3.28b),  one  assumes  that  bn/bt  -  bttj/bt  (i.e.,  the 
net  charge  density  is  approximately  constant  in  time) ,  then  Allis 
and  Rose  (summarized  in  Ref.  13)  have  shown  that  the  effect  of  the 
diffusion  field  E^  on  the  electrons  may  be  accounted  for  by  a  new 
effective  diffusion  constant,  Dfc,  and  the  behavior  of  the  heavy 
ions  then  ignored.  That  is,  they  have  shown  that  the  movement  of 
the  electrons  may  be  described  by  Equation  (3.27a)  but  with 
(~)  •(7>  -  ~ 7  given  by 

n  (c)  -  -  V  (Dtn) 

where  the  new  effective  constant  Dfc  was  obtained  by  numerical 
solution  of  Equations  (3.27a),  (3.27b),  (3.29),  (3.28a),  and  (3.28b) 
with  the  assumption  of  constant  D@,  Dj,  ne,  and  It  is 

believed  that  taking  these  parameters  constant  does  not  greatly 
disturb  the  results. 

The  numerical  results  for  are  given  in  Figure  8  of 

Reference  13  and  may  be  represented  approximately  by  the  formula 


1 


3-23 


l°g(^)-  log(0.9)  l°9<°-9>  [1o9(tt~) 


when  1  <  (D^D^  ■£  (Da/De)  *  Otherwise, 

max  (Dt/De)  ■  1  ,  min  (D^Dg)  -  Da/De 

Here  is  the  Debye  distance  given  by 


log (0.018) 
(3.41a) 

(3.41b) 


V  -  ( 


(3.42) 


D  \  4irne2 

A  is  the  "diffusion  length,"  A  ■  V2rv/n,  and  D=  is  the  ambi- 

CL 

polar  diffusion  constant  given  by 


U_D  +  LL  D_ 
e  ^e  I 


(3.43) 


Equation  (3.41)  states  that  the  logarithm  of  Dt/De  varies 
linearly  from  zero  (free  diffusion  where  Dt  ■  and 
4ttA2/Ad2  ■  0.018)  to  log  Da/De  (full  anibipolar  diffusion  where 
Dt  ■  Da  and  4ttA2/>d2  ■  180) .  Then  three  regions  can  be  distin¬ 
guished 


Free  diffusion: 


Transition: 


Full  ambipolar: 


(' 


j-J  <  1.2X10"3 


Dt  ’  De 


D 


a 


}  (3.44) 
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Since  m  «  Mj,  Equations  (3.33)  and  (3.35)  indicate  |xe  »  p,,. 
and,  using  (3.40),  that  (3.43)  may  be  written  to  good  approximation 

Da  =  (j  2s§^+  kT^  (3.45a) 

With  a  relatively  cold  gas  in  an  applied  field  under  station¬ 
ary  conditions,  the  average  electron  energy  would  be  about  (see, 
e.g.,  Fig.  3.1  and  Ref.  24  for  strong  fields) 

me2/ 2  =  (2/3 )ux  »  kT 

where  ux  is  the  lowest  excitation  potential  and  is  about  an 
electron  volt  for  oxygen.  Then  Equation  (3.45a)  may  be  written 
approximately 


Da  *  (4/9)plUx  (3.45b) 

If  we  take  (3/2)  (mcV2)  55  ux  =  1  ev  in  Equations  (3.42)  and 
(3.44),  we  find  for  the  three  regions 


Free  diffusion:  A2n  £  2.9xl02  cm-1 

Transition  diffusion:  2.9xl02  A2n  3.5x10®  cm' 

Full  ambipolar:  A2n  >  3.5x10°  cm-1 


>  (3.46) 


Thus  it  is  seen  that  in  re-entry  problems  the  diffusion  will  be 
transition  or  full  ambipolar. 

Diffusion  constants  for  ions  and,  therefore,  the  ambipolar 
diffusion  constant  may  be  determined  from  the  polarizability  6 
of  the  neutral  particles  by  use  of  the  formula  (see  Ref.  25, 

Eq.  (14.2  )) 
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3 

kT  (Mj  + 

(£)  <0.0422)r  (3  -i) 

27TMJM 

(3.47) 


where  ■  26(e/c)2,  T  is  the  usual  gamma  function,  6  is  the 

average  polarizability  of  the  neutral  molecules  and  e  is  the 
dielectric  constant. 


3.2.4  The  ionization  oscillation  correction  of  Gould  and  Roberts 
In  Reference  3  high-frequency  breakdown  condition  computations 
were  made  with  the  use  of  d-c  cold-air  data.  The  relationship 
between  the  a-c  and  d-c  parameters  used  was  a  generalization  of 
the  use  of  Equation  (3.19). 

In  Appendix  A,  it  is  noted  that  a  second  harmonic  heating 
contribution  to  the  equation  governing  the  RVDF  is  assumed  negligible 
at  very  high  frequencies.  In  Reference  3  an  attempt  was  made  to 
re-introduce  this  correction  in  a  somewhat  phenomenological  way. 

Gould  and  Roberts  noted  that  if  the  electron  energy  actually  oscil¬ 
lated  (i.e.,  was  modulated)  with  time,  then  some  correction  for 
this  oscillation  should  be  necessary  in  computing  ^v"^,  since 
Vj  depends  strongly  on  the  number  of  electrons  with  energy  above 
the  ionization  potential .  They  computed  the  energy  modulation  by 
use  of  the  macroscopic  energy  equation  with  the  instantaneous  heat¬ 
ing  term  as  given  by  Equation  (3.38)  and  with  energy  loss  per  unit 
time  determined  from  experiment.  With  this  modulated  energy,  the 

-  vR^  is  determined  from 

(V)  5  -  *A  -  ^>-<h(5>  -  V5)]  V5))  (3-48) 


new  average 


-  v, 


where  u  is  the  function  of  time  computed  from  the  energy  equa¬ 
tion,  aT(u)  is  the  d-c  experimentally  determined  ionizations 
per  unit  drift  length  (Townsend  coefficient),  i|a(u)  is  the  d-c 
experimentally  determined  attachments  per  unit  drift  length,  and 
Vj(Ti)  _is  the  d-c  experimentally  determined  drift  velocity, 
=(c),  determined  from  Vd(E/p)  and  "u(E/p). 


L 
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The  (v>  determined  in  this  way  is  plotted  in  Figure  3.2  as 
a  function  of  E7p  for  various  values  of  p X.  The  energy  modu¬ 
lation  correction  goes  to  zero  for  very  high  frequencies  (pk  ■  0) 
and  then  the  relationship  between  a-c  and  d-c  may  be  expressed 
by  Equation  (3.19)  alone.  It  should  be  noted  that  p  used  for 
Figure  3.2  is  not  the  actual  pressure  but  instead 

p*  -  760  -2-  Z  mm  Hg 
Po 

where  Z  is  the  particle  multiplicity  factor  usually  used  in  the 
equation  of  state  in  the  form  p  ■  pZRT,  and  pQ  is  the  sea  level 
density.  In  our  computations  (see  Sections  4  and  7)  p/pQ  is 
about  10-3  to  10-4 ,  Z  -  2,  and  X  =  100  cm.  Then  the  value  of 
p*  to  be  used  is  about  7.2  to  72.  It  is  seen  in  Figure  3.2  that 
the  modulation  correction  here  will  be  considerably  less  than 
1  volt/ mm  Hg  and,  therefore,  no  modulation  correction  is  needed. 

3.3  Data  for  Cold  Air 

In  order  to  check  the  accuracy  of  the  model  used  in  computa¬ 
tion,  it  is  desirable  to  compare  some  computed  results  with  experi¬ 
mental  results  for  the  same  antenna  configuration.  This  is  done 
in  Section  6.  Since  these  experiments  were  conducted  in  cold  air 
under  zero  velocity  conditions,  only  the  high-frequency  ionization 
coefficient  C,  for  cold  air  is  needed. 

3.3.1  The  high-frequency  ionization  coefficient 

In  Reference  5 ,  Figure  2 ,  data  for  v  and  D  have  been 
reduced  and  plotted  in  the  convenient  form  of  £  "  v/DEea  as  a 
function  of  E^p.  The  E/p  axis  of  this  figure  must  be  shifted 
6  volts,  so  that  the  abscissa  reads  about  6  volts  higher,  because 
a  full  energy  modulation  correction  has  been  made  which,  as  pre¬ 
viously  discussed,  does  not  apply  here  (see  Section  3.2.4).  This 
figure  (see  Fig.  3.3)  is  matched  to  within  4-percent  accuracy  by 
the  equation 
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1 


Seoia  ■  <°-243>  (t  ■  30)  xl0'‘  C^tT  (3-49) 

with 


<  110 


volts/ cm-mm  Hg 


This  curve  was  derived  by  expressing  £  in  terms  of  ioniza¬ 
tion  and  attachment  coefficients  as 


3  (a-/ p)  -  (t)  / p) 

2  viyir . 


(3.50) 


The  data  of  (c^/p)  -  (na/p)  were  derived  from  Reference  29a  and 
values  of  u  as  a  function  of  E/p  were  taken  from  References 
29b  and  29c.  Although  some  extrapolation  of  this  data  is  possible 
we  believe  that  our  results  should  not  be  applied  to  situations 
for  which  E^p  is  appreciably  above  110  volts/ cm-mm  Hg. 

3.3.2  The  free  diffusion  coefficient 

For  computing  breakdown  in  gaseous  flows  of  cold  air  one  would 
need  D  in  addition  to  £.  These  cold  air  values  will  also  be 
useful  in  estimating  the  diffusion  for  hot  air.  From  theory,  De 
is  given  by  Equation  (3.40b).  By  using  drift  velocity  measurements 
in  an  applied  field  from  Reference  28,  page  542,  one  can  obtain  |0-e. 
The  average  electron  energy  may  be  taken  from  Reference  2a,  page  83. 
The  obtained  from  these  data  are  plotted  in  Figure  3.4. 

3.4  Data  for  Hot  Air 

Because  of  the  extreme  paucity  of  data  in  hot  air,  the  esti¬ 
mates  here  will  have  to  be  guided  almost  entirely  from  theory. 

For  this  purpose,  the  formulation  for  the  RVDF  and  u  ■  (1/2) me2 
given  by  Allis  (Ref.  13)  will  be  used. 
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3.4.1  Average  energy  and  diffusion  coefficients 

If  in  Equation  (3.23a)  and  (3.23b)  the  diffusion  term  4>  is 
neglected  or  assumed  to  be  about  constant  in  the  elastic  region 
and  if  the  approximation  (3.40a)  is  used  (v0n  -  const),  one  may 
rewrite  the  distribution  function  (3.23)  as 
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(3.51) 


As  pointed  out  in  the  derivation  of  (3.23)  this  equation  does 
not  apply  down  to  c  •  0,  but  all  the  necessary  integrals  in  which 
f  is  used  converge  and  contribute  little  from  the  region  of 
small  c.  From  this  Allis  (Ref.  13)  derives 


where 


where 


J_  u 

1  +  (4wy2l)  +  (2Wq2/63)  •  •  • 

10  uo 

1  +  (w</5) +  Oo*/105)  •  •  • 

(3.52a) 


(3.52b) 


me. 


me. 


and 


Here,  as  mentioned  in  Section  3.1.4,  uo  is  usually  slightly 
larger  than  ux  (the  difference  uQ  -  ux  is  larger  for  larger 
field  strengths) ,  and  then  cQ  is  slightly  larger  than  c*.  The 
term  in  brackets  in  (3.52a)  is  slightly  less  than  unity  and  varies 
very  slowly.  Thus  it  is  seen  that  u  varies  very  slowly  and  will 
be  about  equal  to  or  slightly  above  3u^/10. 
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For  a  choice  of  u^,  it  is  noted  that  the  lowest  excitation 
energy  that  has  been  found  to  be  excited  by  electrons  in  air  or 
air  constituents  was  reported  by  Schulz  in  Reference  30.  Schulz 
reported  a  "...  large  inelastic  process  in  N2  for  electrons  with 
energies  around  2.3  ev  and  that  vibrational  levels  of  the  nitrogen 
molecule,  up  to  v  ■  8,  are  excited."  The  cross  section  for  the 
process  is  not  known. 

If  u  ■  2.3  ev  is  taken  and  if  it  is  assumed  that  W <  1 
and  that  uQ  -  ux  then 

u  a  (2.3)  -  0.7  ev  (3.53) 

Calculations  of  the  electron  energy  distribution  based  upon 
the  2.3  ev  level  reported  by  Schulz  were  made  by  Carleton  and 
McGill  (Ref.  24) .  These  quite  involved  calculations  led  to  mean 
energies  of  0.6  to  1.0  ev  for  fairly  large  values  of  E/p.  This 
agrees  quite  well  with  Equation  (3.53). 

Crompton  and  Sutton  (Ref.-  31)  obtained  average  electron  ener¬ 
gies  in  cold  air  by  measuring  en  E/D  under  conditions  such  that 
mobility  was  balanced  against  diffusion.  Then  from  Equations  (3.28a) 
and  (3.40b)  one  obtains 

u  -  eu(D/ ©M-gE)  -  j  (D/ue) 

Their  results  are  given  in  Figure  3.5.  It  is  seen  that  the  data 
match  a  ux  of  somewhere  between  6  and  8  volts  much  better  than 
0.7  ev  given  by  Equation  (3.53).  It  should  also  be  noted  that 
the  lowest  readily  excitable  (by  electron  impact)  stable  states 
of  the  main  constituents  were  previously  reported  to  be  6.1  ev 
for  and  7.9  ev  for  02  (Ref.  2a).  These  results  imply  that  the 
earlier  mentioned  2.3  ev  is  either  relatively  difficult  to  excite 
or  readily  gives  up  its  energy  in  a  future  impact  or  that  the 
value  of  u  ■  (3/10) ux  cited  by  Allis  is  too  low  for  strong  fields. 
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In  hot  air  still  other  phenomena  may  occur.  Since  the  bound¬ 
ary  layer  may  contain  highly  dissociated  0 2 ,  the  metastable  level 
of  atomic  oxygen  at  1.97  ev  (Ref.  2a)  may  play  an  important  role. 
However,  since  this  possibility  is  uncertain,  we  have  adopted 
ux  ■  6.1  ev  and  used  Equation  (3.52)  to  obtain  Figure  3.6.  Then 
using  Equations  (3.40b)  and  (3.45a),  the  results  shown  in  Figures 
3.7  and  3.8  were  computed  for  the  free  and  ambipolar  diffusion 
coefficients.  The  free-electron  mobilities  for  Equation  (3.40b) 
and  Figure  3.7  were  taken  from  the  collision  frequency  curves  for 
high-temperature  air  in  Reference  15.  The  ion  mobilities  for 
Equation  (3.45a)  and  Figure  3.8  were  computed  from  Equation  (3.47) 
wherein  the  average  molecular  weight  of  the  ions  and  of  the  neu¬ 
trals  were  taken  from  Reference  33  and  the  average  polarizabilities 
for  hot  air  were  taken  from  Reference  34. 

3.4.2  The  effective  ionization  rate 

During  preparation  of  this  report  we  have  found  virtually  no 
data  for  Vj  or  vA  in  hot  air.  Presumably  the  number  of  parti¬ 
cles  making  their  way  to  energies  that  are  sufficiently  high  to 
cause  ionization  would  be  increased  by  the  higher  temperatures 
because  of  the  decreased  ability  to  transfer  this  energy  to  the 
heavy  particles.  However,  this  effect  should  be  slight  because 
the  largest  part  of  the  energy  is  lost  in  inelastic  collisions  and 
almost  all  particles  will  be  in  their  ground  states. 

Another  possibility  at  temperatures  of  2000°  to  4000°  K  and 
pressures  of  around  1  mm  Hg  is  that  the  presence  of  NO  in  significant 
amounts  may  decrease  the  effective  energy  necessary  to  cause  ioniza¬ 
tion  because  of  its  lower  ionization  potential.  However,  breakdown 
measurements  in  relatively  pure  NO  (Ref.  2b)  indicate,  surprisingly, 
that  it  is  actually  slightly  more  difficult  to  break  down  NO  than 
pure  Nj  or  pure  O2 .  Therefore  we  have  assumed  Vj  -  vR  -  vA  to  be 
unchanged  in  hot  air  relative  to  its  cold  air  value.  This  involves 
using  the  cold  air  result  of  Equation  (3.49)  to  obtain 
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(3.54) 


3.5  Mixtures,  Impurities,  and  Use  of  Experimental  Data 

The  effects  of  small  amounts  of  impurities  or  admixtures  in 
certain  gases  have  been  discussed  by  Brown  (Ref.  1)  and  he  shows 
that,  for  selected  gases,  small  concentrations  of  gases  added  to 
the  main  gas  can  cause  large  changes  in  breakdown  conditions.  The 
effect  of  the  impurity  gas  has  not  always  been  well  understood  but 
the  following  phenomena  are  known  to  occur: 

(a)  The  impurity  may  have  a  low,  easily  excited  excitation 
energy  and  will  thus  strongly  influence  the  electron  RVDF.  The 
relationship  of  the  RVDF  to  pertinent  macroscopic  properties  was 
discussed  in  Sections  3.1.4,  3.2.2,  and  3.4. 

(b)  The  presence  of  a  low  ionization  potential  species  in 
the  mixture  may  greatly  increase  the  ionization  coefficient.  This 
was  why,  in  Section  3.4.2,  the  presence  of  NO  in  air  was  expected 
to  increase  the  ionization  rate.  The  fact  that  pure  NO  was  more 
difficult  to  break  down  than  02  and  N2  is  probably  due  to  the  phe¬ 
nomenon  (a)  above.  Since  NO  has  a  strong  permanent  dipole  moment, 
it  is  probably  much  easier  to  excite  its  low  states  by  electron 
impact  than  it  would  be  to  excite  02  or  N2 . 

(c)  High-energy  metastable  states  are  more  difficult  to 
excite  in  the  presence  of  small  quantities  of  low-ionization  energy 
impurities  such  as  Hg. 

(d)  Ramsauer  effects  in  the  impurity  gas  or  changes  in  the 
velocity  distribution  (as  in  (a)  above)  induced  by  the  impurity 
gas  may  have  large  influence  upon  the  electron  mobility  and  dif¬ 
fusion  coefficient.  As  a  possible  example  of  this  phenomenon, 
consider  the  effect  upon  mobility  of  adding  small  amounts  of 
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nitrogen  to  argon  compared  to  the  pure  electron  mobilities  as 
depicted  in  Reference  2a,  Figures  3.8,  3.9,  and  3.12,  when  addi¬ 
tion  of  1  percent  of  N2  to  A  results  in  a  sixfold  increase  in 
mobility  over  pure  argon  at  E/p  ■  1  volt/cm-mm  Hg.  The  mobility 
is  also  more  than  twice  that  of  pure  nitrogen. 

(e)  Attachments  and  recombinations  may  be  strongly  affected 
by  the  presence  of  certain  species. 

Thus,  it  is  seen  that  macroscopic  properties  pertinent  in 
breakdown  may  be  strongly  altered  by  the  presence  of  a  small 
amount  of  impurity  or  by  the  presence,  in  sufficient  amounts,  of 
a  specie  formed  only  at  high  temperatures.  Because  of  these 
effects,  measurements  of  breakdown  in  impure  hot  gases  produced 
in  devices  such  as  arc  jets  may  have  little  meaning. 
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igure  3.1.-  Schematic  of  the  shape  of  the  random  velocity  distribute 
function  for  three  different  field  strengths.  This  distribution  fun< 
tion  is  normalized  to  one  when  multiplied  by  4trc£dc  and  integrated 
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Figure  3.2.-  The  net  ionization  frequency  (v  (E^p)) 
for  various  wavelengths  A.  Taken  from  Reference  3. 


3-36 


o 

■ 

■ 

■ 

■ 

H 

Electron  energy 
~  extrapolated 

■ 

5  60  65  70  75  80  85  90  95  10 

i-mm  Hg 

■ 

■ 

■ 

■ 

H 

■ 

■ 

a 

a 

a 

■ 

a 

■ 

■ 

S 

H 

■ 

■ 

a 

■ 

H 

(Ref.  28,  p.  542) 
r  (Ref.  2a,  p.  83) 

■ 

■ 

■ 

■ 

■ 

■ 

i 

■ 

■ 

■ 

0  5  10  15  20  25  30  35  40  45  50  5 

Eg/p,  volts/cn 

B 

■ 

■ 

g 

■ 

■ 

0. 

*  s 

•n  « 

B 

H 

■ 

■ 

w  y  — 

=L  0  4) 

Nh  o>  c 

>  0 
■  u 

■p  p 
ft  >P  0 

■ 

■ 

■ 

n 

■ 

> 

•rt 

M 

a 

w 

H 

w 

B 

■ 

■ 

■ 

ii 

■ 

■ 

■ 

■  • 

n 

4) 

a 

■ 

■ 

■ 

■ 

ii 

■ 

■ 

■ 

■ 

o 

u 

3 

0 

a 

a 

Data  s 

■ 

■ 

m 

■ 

■ 

■ 

■ 

a 

IS 

— 

— . 

a 

a 

a 

OGOvO^rtO  COV0*NOCO 
•  •••••••••  •• 

s_OTX  &H  uiux-o©9/3uio  ‘d^a 


Figure  3.4.-  Free  electron  diffusion  coefficient  for  cold  air 
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Figure  3.5.-  Average  energy  for  electrons  in  cold  air  versus 
(Eg/p) 2 ,  data  from  Crompton  and  Sutton  (Ref.  31). 
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Figure  3.6.-  Average  energy  of  electrons  as  a  function  of  E  for 
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Figure  3.6.-  Concluded 


3-42 


EBSP 

ft 

6000°  K 

p/p0  -  16“*  . 

r 

\i 

m 

■ 

■■■ 

mm 

m 

t 

i 

• i 

1 

800( 

3°.JS 

p/p 

»■* 

r2 

1 

6000°  K 

p/po  -  10- 

2000  K 

p/p0  -  - 

if? 

■M 

a  / 

/ 

/ 

/ 

i 

i 

i 

E 

e  -  E 

rms 

v2 

f 

2 

v2  + 

jo2/ 

L 

0  10  20  30  40  50  60  70  80  90  100  110  120  130 


Effective  field  E_,  volts/ cm 

© 

Figure  3.8.-  Ambipolar  diffusion  coefficient  in  air 
as  a  function  of  Eg. 
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4.  PROPERTIES  OF  TYPICAL  AERODYNAMIC  FLOW  FIELDS 

The  local  dynamical  and  thermodynamical  properties  of  typical 
aerodynamic  flow  fields  surrounding  reentry  bodies  are  required  as 
inputs  for  the  breakdown  computations.  In  this  section,  we  present 
the  results  of  computations  for  two  simplified  geometrical  shapes 
involving  a  wedge  and  a  sharp-nosed  cone.  Both  the  inviscid  and 
viscous  regions  of  the  flow  field  have  been  treated  including  both 
laminar  and  turbulent  regions.  Assumptions  made  involve:  (1)  The 
use  of  equilibrium,  thermodynamic  properties  for  the  air  in  the 
flow,  (2)  neglect  of  foreign  species  induced  by  ablation  or  trans¬ 
piration  into  the  flow,  (3)  neglect  of  displacement  of  effective 
body  surface  location  by  the  boundary  layer  in  the  inviscid  region, 
and  (4)  negligible  shock  front  thickness. 

Quantities  calculated  for  these  regions  include  local  velocity, 
temperature,  density,  enthalpy  for  both  the  inviscid  and  viscous 
regions  as  well  as  electron  concentration,  weak  field  electron 
collision  frequency  and  strong  field  electron  collision  frequency 
for  the  viscous  boundary  layer.  Electron  concentrations  and  col¬ 
lision  frequencies  are  computed  from  the  thermodynamic  state  of 
the  gas.  These  calculations  have  been  incorporated  into  the  com¬ 
puter  program  for  the  viscous  flow  field,  as  described  in  Reference  36. 


4.1  Inviscid  Flow-Field  Calculations 

4.1.1  Inviscid  flow  field  -  wedges 

The  inviscid  flow  fields  surrounding  a  two-dimensional  wedge 
have  been  calculated  using  the  method  of  W.  E.  Moeckel  (Ref.  37). 
The  equations  used  in  these  calculations  are: 

P  p  v  2 

~  “  (1  -  u  )  ~l  —  sin2  8+1  (4.1) 

P,  r  p 


(4.2) 
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tan  (0  -  A)  ■  ur  tan  0  (4.3) 

where  subscript  l's  refer  to  free-stream  quantities  and  sub¬ 
script  2's,  to  quantities  in  the  flow  field  and 


is  the  ratio  of  the  components  of  velocities  and  ,  normal 

to  the  oblique  shock  wave.  The  quantity,  A  refers  to  the  angle 
between  the  surface  of  the  wedge  and  the  direction  of  the  free- 
stream  velocity  vector,  and  6  is  the  oblique  shock  angle  with 
respect  to  the  direction  of  the  same  vector. 

The  equations  above  are  solved  by  iteration,  using  various 
input  conditions.  The  ARDC  19b9  Atmosphere  (Ref.  38)  has  been 
used  throughout  this  work,  as  well  as  for  the  conical  cases. 

The  results  of  the  calculations  for  the  inviscid  flow  fields 
are  presented  in  Tables  4.1(a)  through  (d) .  Because  stream  lines 
in  the  flow  field  are  parallel  to  the  body  surface,  the  quantities 
plotted  are  constant  throughout  the  flow  field  and  so  represent 
the  values  next  to  the  body  surface. 

4.1.2  Inviscid  flow  field  -  sharp-nosed  cones 

The  inviscid  flow-field  parameters  associated  with  the  hyper¬ 
sonic  flight  of  a  sharp-nosed  cone  at  zero  angle  of  attack  have 
been  derived  from  the  work  of  Romig  (Ref.  39)  using  the  ARDC  1959 
atmospheric  data  (Ref.  38).  Her  work  assumes  continuum  flow, 
negligible  boundary-layer  thickness,  negligible  shock-front  thick¬ 
ness,  and  steady-state  flow.  The  results  of  these  calculations 
for  a  number  of  typical  cases  are  given  in  Tables  4.2(a)  through  (d) . 
The  quantities  plotted  are  those  calculated  for  the  region  adjacent 
to  the  body  surface. 
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4.2  Boundary-Layer  Flow-Field  Calculations 

The  viscous  flows  associated  with  each  of  the  geometric  shapes 
described  in  the  previous  sections  have  been  calculated  using  the 
method  of  Reference  36.  This  reference  describes  a  computational 
procedure  and  computer  program  based  upon  a  method  proposed  by 
Rubesin  (Ref.  34).  A  detailed  discussion  of  the  procedure  is 
given  in  Reference  36.  The  inviscid  flow  fields  discussed  in  the 
previous  sections  furnish  the  inputs  to  this  program. 

In  the  case  of  cones,  the  boundary -layer  properties  were  cal¬ 
culated  at  distances  back  from  the  apex  along  the  surface  such 
that  the  antennas  used  would  fit  into  the  space  available.  These 
distances  along  the  surface  occurred  at  distances  of  5.73,  3.82, 
and  2.86  feet  for  the  10°,  15°,  and  20°  cones,  respectively. 

Sample  results  of  the  application  of  this  program  to  the 
flow-field  calculation  are  given  in  Figures  4.1  through  4.7. 
Complete  results  of  these  calculations  are  available  from  Vidya, 
Inc.  However,  since  the  results  of  these  calculations  are  quite 
lengthy  for  the  number  of  cases  and  the  number  of  variables  con¬ 
sidered,  only  sample  results  are  presented.  Furthermore,  in  apply¬ 
ing  the  results  of  these  calculations  to  the  breakdown  calculation, 
certain  approximations  have  been  used,  as  are  discussed  in  the 
next  section. 

4.3  Representation  of  the  Boundary-Layer  Profile  with  the 

Function  B 

The  calculations  of  the  previous  section  are  now  used  to 
obtain  values  for  the  quantity  m  \f/D,  defined  in  Section  2.1.2. 
Two  formulations  for  breakdown  computations  discussed  in  the  next 
section  involve  the  use  of  (a)  a  mean  constant  value  for  the  dif¬ 
fusion  coefficient  D,  and  (b)  a  variable  diffusion  coefficient 
which  is  a  function  of  the  local  conditions  at  each  point  through¬ 
out  the  flow  field.  As  described  in  Section  7,  the  computations 
of  breakdown  have  been  done  initially  for  the  constant  D  case. 
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The  velocity  profile  vf(y)  was  approximated  by  the  function 

V  -  (1  -  e“by)  (4.5) 

where  the  parameter  b  was  adjusted  to  give  the  best  fit  to  the 
velocity  profile  in  the  boundary  layer,  and  is  given  by  the 

velocity  at  the  edge  of  the  boundary  layer.  These  data  are  tabu¬ 
lated  in  Table  4.3. 

The  values  for  D  were  taken  at  a  point  in  the  boundary 
layer  corresponding  to  the  position  at  which  the  electric  field 
had  fallen  to  one-half  its  value  at  z  -  0  on  the  surface  (see 
also  Section  7).  The  antenna  for  which  the  calculations  have  been 
performed  is  a  slot  1/2  inch  wide  by  17  inches  long,  oriented  on 
the  surface  with  its  width  normal  to  the  gaseous  flow  direction. 

The  near  field  of  this  antenna  has  been  examined  by  Scharfman  and 
Morita  (Ref.  4).  As  discussed  in  Section  6.1,  their  results  show 
that  the  field  at  z  »  0  may  be  approximated  by  E  *  Eq  e-€iy 
where  e l  ■  0.540  cm-1,  so  that  at  E/Eq  *  1/2, 
y  ■  1.293  cm  «  0.0425  feet. 

The  values  for  the  temperature,  density,  and  electron  density 
on  the  boundary  layer  were  obtained  at  y  -  0.0425  feet,  and  values 
of  D  were  obtained  from  Figures  3.7  and  3.8.  These  values  are 
listed  in  Table  4.3  for  each  of  the  flow  situations  considered. 
Finally,  values  of  Bq  -  VQ/D  are  also  given  in  Table  4.3. 
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TABLE  4.1.-  INVISCID  FLOW-FIELD  CALCULATIONS  FOR  CONES. 

(a)  h  -  80,000  ft. 

T  -  390°  R,  p  -  2.79X10"3  lb/ft3,  pi  -  58.1  lb/ft2,  hx  -  93.9  BtVlb 
hwall  "  646  Btu/lb»  Twan  “  2500°  K,  Prandtl  No.  -  0.75 


A 

V 

V 

p 

p 

h 

T 

l 

2 

'2 

2 

2 

2 

(deg) 

(ft/sec 

xlO-3) 

(ft/aec 

XlO"3) 

(lb/ft3 

XlO2) 

(lb/ ft2) 

(Bti/lb) 

(°R) 

10 

25 

24.393 

1.537 

1820 

530 

2028 

15 

25 

23.77 

1.942 

4011 

1020 

3685 

20 

25 

22.804 

2.417 

6684 

1440 

4875 

10 

15 

- 

1.145 

726.6 

235 

1014 

15 

15 

14.142 

1.425 

1511 

415 

1677 

20 

15 

13.675 

1.676 

2616 

678 

2613 

10 

10 

- 

.950 

372 

150 

663 

15 

10 

- 

1.145 

726.6 

320 

975 

20 

10 

8.832 

1.327 

1162 

345 

1404 

I 


TABLE  4.1.-  CONTINUED, 
(b)  h  -  120,000  ft. 


Tx  -  451.37°  R,  pt  -  4.0851X10"4  lb/ft3,  p  -  9.8372  lb/ft2 
hi  -  1. 085X1 O2  Btu/lb,  bwall  -  646  Btu/lb,  T  -  2500°  K 

Prandtl  No.  ■  0.75 


A 

V 

V 

p 

p 

h 

T 

1 

2 

r2 

2 

2 

2 

(deg) 

(ft/sec 

xlO'3) 

(ft/sec 

xlO-3) 

(lb/ft3 

XlO3) 

(lb/ ft2) 

(Btu/lb) 

(°R) 

10 

35 

34.321 

2.63 

501.7 

930 

3430 

15 

35 

33.541 

3.74 

1151 

1710 

5281 

20 

35 

32.434 

4.40 

1869 

3165 

6996 

10 

25 

24.454 

2.14 

275.4 

625 

2076 

15 

25 

23.979 

2.74 

570.6 

1055 

3792 

20 

25 

23.238 

3.53 

1003 

1600 

5055 

10 

15 

- 

1.70 

102.3 

240 

1038 

15 

15 

14.177 

2.02 

216.4 

430 

1738 

20 

15 

13.820 

2.41 

373.8 

690 

2686 
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1 

I 

I  TABLE  4.1.-  CONTINUED. 

I  j  (c)  h  -  200,000  ft. 

1 


T  -  449°  R,  p  -  1.97X10”5  lb/ft3,  px  -  0.472  lb/ft2  ,  -  108  Btu/lb 

hwall  "  646  Btu/lb»  Twall  “  2500°  K»  Prandtl  No.  -  0.75 


A 

V 

V 

p 

p 

h 

T 

l 

2 

*2 

2 

2 

2 

(deg) 

(ft/sec 

xlO-3) 

(ft/ sec 
xlO-3) 

(lb/ft3 

XlO4) 

(lb/ ft2) 

(BtVlb) 

(°R) 

10 

45 

44.16 

1.6 

40.08 

1210 

4131 

15 

45 

43.19 

2.2 

84.87 

2925 

5927 

20 

45 

41.71 

2.4 

150.9 

4590 

8531 

10 

35 

34.35 

1.3 

24.05 

1080 

3772 

15 

35 

33.54 

1.9 

55.17 

1590 

4714 

20 

35 

32.43 

2.2 

91.94 

3130 

6151 

10 

25 

24.45 

1.1 

13.20 

525 

2065 

15 

25 

23.94 

1.4 

27.35 

1080 

3772 

20 

25 

23.13 

1.8 

48.57 

1500 

4580 

10 

10 

- 

.7 

2.829 

190 

727.4 

15 

10 

- 

.8 

5.187 

230 

1033 

20 

10 

- 

1.0 

8.393 

360 

1459 
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TABLE  4.1.-  CONCLUDED, 
(d)  h  -  250,000  ft. 


T  -  328.2°  R,  p  -  2.493x10"®  lb/ ft3,  -  0.0436  lb/ ft2 
h  -  78.9  Btv/lb,  hwall  -  646  Bti^lb,  T%(all  -  2500°  K 

Prandtl  No.  ■  0.75 


A 

V 

V 

p 

p 

h 

T 

i 

2 

'£ 

2 

2 

2 

(deg) 

(ft/  sec 

X10"3) 

(ft/ sec 
X10"3) 

(lb/ft3 

xlO5) 

(lb/ ft2) 

(Btu/lb) 

(°R) 

10 

15 

- 

1.1 

0.6110 

230 

919 

15 

15 

14.07 

1.4 

1.309 

387 

1575 

20 

15 

13.42 

1.6 

2.182 

685 

2461 

10 

25 

24.31 

1.5 

1.658 

490 

1969 

15 

25 

23.56 

2.2 

3.404 

750 

2855 

20 

25 

22.87 

2.9 

6.022 

950 

3446 

10 

35 

34.10 

2.0 

3.011 

720 

2757 

15 

35 

33.02 

3.0 

6.546 

1160 

3545 

20 

35 

31.46 

3.1 

11.35 

2760 

5382 
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TABLE  4.2.-  INVISCID  PLCW-PIELD  CALCULATIONS  POR  WEDGES. 

(a)  h  -  80,000  ft. 


T  -  390°  R,  p  -  2.79X10"3  lb/ ft3,  pi  -  58.1  lb/ ft2,  -  93.9  BtVlb 
hwall  "  646  Btu/lb»  Twaii  "  2500°  K»  Prandtl  No.  -  0.75 


A 

V 

e 

V 

p 

p 

h 

T 

M 

l 

2 

r2 

*2 

2 

2 

2 

(deg) 

(ft/ sec 
xlO-3) 

(deg) 

(ft/sec 

xlO-3) 

(lb/ ft3 
XlO2) 

(lb/ ft2) 

(Btu/lb) 

(°R) 

Mach 

No. 

10 

25 

12.2 

24.45 

1.57 

2052 

635 

2450 

10.4 

15 

25 

17.8 

23.8 

1.84 

4354 

1280 

4420 

B ■ 

20 

25 

23.1 

23.0 

2.22 

7343 

2030 

5987 

10 

15 

13.0 

14.63 

1.22 

824 

375 

1266 

8.5 

15 

15 

18.6 

14.24 

1.50 

1674 

535 

2102 

6.52 

20 

15 

24.3 

13.71 

1.68 

2814 

690 

3136 

5.18 

10 

10 

14.3 

9.72 

.953 

431 

210 

848 

6.84 

15 

10 

19.7 

9.45 

1.22 

818 

310 

1260 

5.50 

20 

10 

25.4 

9.08 

1.41 

1336 

435 

1778 

4.49 

TABLE  4.2.-  CONTINUED, 
(b)  h  -  120,000  ft. 


T  -  451.37°  R,  pi  -  4.0851X10-4  lb/ft3,  pi  -  9.8372  lb/ft2 
hi  -  1 . 085x1 02  Btu/lb,  hwall  -  646  Btu/lb,  T  n  -  2500°  K 

Prandtl  No.  -  0.75 


A 

V 

e 

V 

p 

p 

h 

T 

M 

i 

2 

r2 

2 

2 

2 

2 

(deg) 

(ft/sec 

XlO"3) 

(deg) 

(ft/ sec 
XlO"3) 

(lb/ ft3 
XlO2) 

(lb/ ft2) 

(Btu/lb) 

(°R) 

Mach 

10 

35 

11.9 

34.3 

2.61 

566 

1150 

4056 

11.6 

15 

35 

17.0 

33.4 

3.57 

1188 

2200 

5853 

9.5 

20 

35 

22.2 

32.4 

4.26 

2026 

3520 

7544 

7.49 

10 

25 

12.3 

24.4 

2.26 

3  02 

645 

2514 

10.3 

15 

25 

17.8 

23.8 

2.71 

637 

1360 

4392 

7.8 

20 

25 

22.9 

23.1 

3.39 

1068 

2020 

5620 

6.72 

10 

15 

13.2 

14.6 

1.73 

123 

372 

1337 

8.28 

15 

15 

18.7 

14.2 

2.14 

247 

550 

2167 

6.42 

20 

15 

24.4 

13.7 

2.43 

414 

850 

3197 

5.13 
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TABLE  4.2.-  CONTINUED, 
(c)  h  -  200,000  ft. 


T  -  449°  R,  px  -  1.97X10"5  lb/ ft3,  pi  -  0.472  lb/ ft2,  hx  -  108  BtVlb 
hwall  "  646  Btu/lb*  Twall  "  2500°  K»  Pr»ndtl  No.  -  0.75 


A 

V 

6 

V 

^2 

p 

h 

T 

M 

l 

2 

*2 

2 

2 

2 

(deg) 

(ft/ sec 
XlO"3) 

(deg) 

(ft/sec 

XlCT3) 

(lb/ft3 
XI O4) 

(lb/ ft2) 

(Btu/lb) 

(°R) 

Mach 

10 

45 

11.4 

44.1 

1.62 

43.0 

- 

4770 

14.1 

15 

45 

16.5 

43.2 

2.20 

91.7 

- 

6568 

10.6 

20 

45 

21.8 

41.8 

2.45 

158.0 

- 

9365 

8.63 

10 

35 

11.8 

34.27 

1.29 

27.2 

1130 

3938 

11.9 

15 

35 

16.8 

33.52 

1.87 

56.6 

1383 

5204 

10.1 

20 

35 

22.1 

32.4 

2.18 

97.0 

2870 

6974 

7.64 

10 

25 

12.2 

24.4 

1.09 

14.5 

646 

2510 

10.3 

15 

25 

17.6 

23.9 

1.35 

30.5 

1240 

4179 

8.14 

20 

25 

22.6 

23.1 

1.77 

50.8 

1950 

5032 

7.13 

10 

10 

14.5 

9.71 

.644 

3.15 

220 

916 

6.58 

15 

10 

19.9 

9.44 

.829 

5.88 

320 

1331 

5.36 

20 

10 

25.5 

9.07 

.970 

9.53 

452 

1844 

4.41 
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TABLE  4.2.-  CONCLUDED, 
(d)  h  -  250,000  ft. 


T  -  328.2°  R,  pA  -  2.493X10-6  lb/ ft3,  Pa  -  0.0436  lb/ ft2 
h  -  78.9  Btu/lb.  bwall  -  646  Btu/lb,  Twall  -  2500°  K 

Prandtl  No.  ■  0.75 


A 

V 

e 

V 

p 

p 

h 

T 

M 

l 

2 

*2 

2 

2 

2 

2 

(deg) 

(ft/ sec 
xlO-3) 

(deg) 

(ft/sec 

XlO*3) 

(lb/ft3 

XlO6) 

(lb/ ft2) 

(BtVlb) 

(°R) 

Mach 

NO. 

10 

15 

12.9 

14.640 

1.130 

0.7205 

293 

1197 

8.73 

15 

15 

18.5 

14.251 

1.363 

1.480 

515 

2037 

6.62 

20 

15 

24.2 

13.718 

1.527 

2.496 

820 

3066 

5.26 

10 

25 

12.2 

24.456 

1.427 

1.817 

615 

2389 

10.5 

15 

25 

17.5 

23.863 

1.786 

3.821 

1203 

3955 

8.44 

20 

25 

22.5 

23.124 

2.385 

6.380 

1875 

4654 

7.43 

10 

35 

11.76 

34.281 

1.689 

3.405 

1095 

3755 

12.34 

15 

35 

16.69 

33.541 

2.529 

7.102 

2045 

4798 

10.52 

20 

35 

23.03 

32.464 

2.846 

12.23 

3580 

6703 

7.81 
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TABLE  4.3.-  INPUT  DATA  FOR  ANTENNA  BREAKDOWN  COMPUTATIONS. 

(a)  Cones. 


Case 

N.  No. 

Parameters  n. 

2 

3 

4 

5 

6 

Altitude 

(ft) 

250,000 

250,000 

250,000 

250,000 

250,000 

Velocity 
(ft/  sec) 

15,000 

15,000 

25,000 

25,000 

25,000 

Body  half  angle 
(deg) 

15 

20 

10 

15 

20 

VQ  (cn/aec) 

4.58X105 

4.58x10s 

7.63X105 

7.63x10s 

7.63x10s 

D  (cm2/sec) 

1.1x10s 

5.8X104 

1.5x10s 

6.5X104 

2.4X104 

b  (cm-1 ) 

0.56 

0.699 

0.495 

0.784 

1.21 

BQ  (cm-1) 

4.16 

7.90 

5.09 

11.7 

31.8 

Z 

1.10 

1.10 

1.21 

1.20 

1.17 

Temperature 

2320 

2340 

3550 

3280 

2665 

Pressure 
(mm  Hg) 

0.4695 

0.785 

0.595 

1.220 

2.160 

Density 

<P/PC> 

6.88X10-5 

1.17X10-4 

4.98X10-5 

1.13X10-4 

2. 79x1 0-4 

VDa 

89 

90 

87.5 

84 

90 
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TABLE  4.3.-  CONTINUED, 
(a)  Continued. 


Case 

N.  No. 

7 

8 

9 

10L 

10T 

Parameters'^. 

Altitude 

(ft) 

250,000 

250,000 

250,000 

200,000 

200,000 

Velocity 
(ft/  sec) 

35,000 

35,000 

35,000 

45,000 

45,000 

Body  half  angle 
(deg) 

10 

15 

20 

10 

10 

VQ  (on/ sec) 

1.068X10® 

1.068x10® 

1.068X10® 

1.37X10® 

1.37X10® 

D  (cm2/ sec) 

1.3X105 

5.5X104 

3 . OxlO4 

3.7X103 

1 .OxlO4 

b  (cm-1 ) 

0.620 

1.022 

1.288 

2.93 

3.87 

Bq  (cm’1) 

8.23 

19.45 

35.6 

371 

137 

Z 

1.35 

1.32 

1.24 

1.04 

1.53 

Temperature 

(°K) 

4660 

4555 

4135 

2300 

5120 

Pressure 
(mm  Hg) 

1.081 

2.345 

4.075 

14.38 

14.38 

Density 

<P/PG> 

6.12X10"5 

1.46X10’4 

2.87X10’4 

2.22X10’3 

8.29X10’4 

VD, 

72.5 

77 

87 

100 

84.5 
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TABLE  4.3.-  CONTINUED, 
(a)  Continued. 


N.  Case 

N.  No. 

11L 

11T 

12L 

12T 

13L 

Parameters^''*.. 

Altitude 

(ft) 

200,000 

200,000 

200,000 

200,000 

200,000 

Velocity 
(ft/ sec) 

45,000 

45,000 

45,000 

45,000 

35,000 

Body  half  angle 
(deg) 

15 

15 

20 

20 

10 

VQ  (cn/*«c) 

1.37X10® 

1.37X10® 

1.37x10® 

1.37X10® 

1.068X10® 

D  (cm2/sec) 

3 .lxlO3 

4.6X103 

2.5X103 

2.5X10® 

5.5X10® 

b  (cm-1) 

4.03 

4.14 

4.77 

4.53 

1.86 

Bq  (cm-1) 

442 

298 

548 

548 

194.5 

Z 

1.12 

1.24 

1.25 

1.25 

1.05 

Temperature 

(°K) 

3290 

4720 

4740 

4740 

2090 

Pressure 
(mm  Hg) 

30.4 

30.4 

54.1 

54.1 

8.65 

Density 

(P/P0> 

2.87X10-3 

1.87X10-3 

3.03X10-® 

3.03X10-® 

1. 45X10-® 

VD, 

94 

83 

84 

94 

99.5 
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TABLE  4.3.-  CONTINUED, 
(a)  Continued. 


n.  Case 

N.  No. 

13T 

14L 

14T 

15L 

15T 

Parameters 

Altitude 

(ft) 

200,000 

200,000 

200,000 

200,000 

200,000 

Velocity 

(ft/sec) 

35,000 

35,000 

35,000 

35,000 

35,000 

Body  half  angle 
(deg) 

10 

15 

15 

20 

20 

VQ  (cnv/sec) 

1.068X10® 

1.068X10® 

1.068x10® 

1.068X10® 

1.068x10® 

D  (cm2/s®c) 

l.lXlO4 

2.2X103 

3.2X103 

2.7X103 

2.7X103 

b  (cm-1) 

3.9 

3.28 

6.95 

3.88 

7.48 

Bq  (cm-1) 

97 

486 

334 

396 

396 

Z 

1.21 

1.06 

1.10 

1.17 

1.17 

Temperature 

(°K) 

3580 

2620 

2980 

3417 

3417 

Pressure 
(mm  Hg) 

8.65 

19.8 

19.8 

32.9 

32.9 

Density 

(P/P0} 

7.35X10-4 

3.80X10-3 

2.18X10-3 

2.94X10-3 

2.94X10-3 

VDa 

87 

99 

95.5 

94 

94 
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TABLE  4.3.-  CONTINUED, 
(a)  Concluded. 
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TABLE  4.3.-  CONTINUED, 
(b)  Wedges. 


Nv  Case 

No. 

19 

20 

21 

22 

23 

Parameters^-^^ 

Altitude 

(ft) 

250,000 

250,000 

250,000 

250,000 

250,000 

Velocity 
(ft/ sec) 

15,000 

15,000 

15,000 

25,000 

25,000 

Body  half  angle 
(deg) 

10 

15 

20 

10 

15 

Vq  (cn/sec) 

4.58X105 

4.58x10s 

4.58x10s 

7.63x10s 

7.63x10s 

D  (cm2/ sec) 

1.8X105 

8Xl04 

5xl04 

9.5X104 

7.9X10“ 

b  (cm-1 ) 

0.244 

0.341 

0.462 

0.341 

0.525 

Bq  (cm-1) 

2.54 

5.7 

9.2 

8.04 

9.67 

Z 

1.05 

1.09 

1.06 

1.21 

1.22 

Temperature 

(°K) 

2167 

2310 

2427 

3333 

3890 

Pressure 
(mm  Hg) 

0.261 

0.531 

0.895 

0.652 

1.37 

Density 

<P/  PG) 

4.17X10"5 

7 . 89xl0-s 

1.26X10"4 

5.85xl0“s 

1.06x10““ 

VDa 

88 

89 

89 

- 

82 

79.5 

TABLE  4.3.-  CONTINUED 


(b)  Continued. 


v —  Case 

Xv  No. 

24 

25 

26 

27 

31L 

Parameters^^^ 

Altitude 

(ft) 

250,000 

250,000 

250,000 

250,000 

200,000 

Velocity 

(ft/sec) 

25,000 

35,000 

35,000 

35,000 

35,000 

Body  half  angle 
(deg) 

20 

10 

15 

20 

10 

VQ  (cn/aec) 

7.63X105 

1.068X10® 

1.068X10® 

1.068X10® 

1.068X10® 

D  (cm2/sec) 

4.5X104 

1.2x10s 

5.6X104 

4.3X104 

4.7X103 

b  (cm-1 ) 

0.644 

0.414 

0.614 

0.772 

2.01 

Bq  (cm-1) 

17 

8.91 

19.1 

24.8 

227 

Z 

1.23 

1.34 

1.40 

1.42 

1.02 

Temperature 

(°K) 

4060 

4580 

4830 

5060 

2130 

Pressure 
(mm  Hg) 

2.29 

1.22 

2.55 

4.39 

9.75 

Density 

(p/Po* 

1.68X10"4 

6.97xl0“5 

1.35xl0“4 

2.24X10"4 

1.58X10“® 

VD. 

80 

73 

75 

75 

100 
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TABLE  4.3.-  CONTINUED, 
(b)  Continued. 


N.  Case 

n.  No. 

31T 

32L 

32T 

33L 

33T 

Parametersv^N^ 

Altitude 

(ft) 

200,000 

200,000 

200,000 

200,000 

200,000 

Velocity 

(ft/sec) 

35,000 

35,000 

35,000 

35,000 

35,000 

Body  half  angle 
(deg) 

10 

15 

15 

20 

20 

VQ  (cnv/sec) 

1.068X10® 

1.068X10® 

1.068x10® 

1.068x10® 

1.068X10® 

D  (cm2/ sec) 

1.15X104 

2.6X103 

6X10  3 

2.5X103 

4x10  3 

b  (cm-1 ) 

1.57 

3.17 

2.12 

3.47 

2.13 

Bq 

93 

414 

178 

427 

267 

Z 

1.27 

1.07 

1.23 

1.18 

1.25 

Temperature 

(°K) 

4670 

2890 

4500 

3870 

4830 

Pressure 
(mm  Hg) 

9.75 

20.3 

20.3 

34.8 

34.8 

Density 

<P/Po) 

5.76X10"4 

2.90X10"3 

1.33X10-3 

3.31X10-3 

1.95X10-3 

VDa 

79.5 

97 

83 

91 

82 
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TABLE  4.3.-  CONTINUED, 
(b)  Continued. 


NV  Case 

No. 

34L 

34T 

35L 

35T 

36L 

Parameters^. 

Altitude 

(ft) 

200,000 

200,000 

200,000 

200 , 000 

200,000 

Velocity 

(ft/sec) 

25,000 

25,000 

25,000 

25,000 

25,000 

Body  half  angle 
(deg) 

10 

10 

15 

15 

20 

VQ  (cn\/sec) 

7.63X105 

7.63X105 

7.63x10s 

7.63x10s 

7.63x10s 

D  (cm2/aec) 

7.8X103 

1.35X104 

4.5x10s 

6.8x10s 

3.8X103 

b  (cm-1) 

1.86 

1.78 

2.09 

3.59 

2.57 

Bq  (cm"1) 

98 

56.6 

171 

112 

206 

Z 

1.00 

1.15 

1.07 

1.13 

1.07 

Temperature 

(°K) 

1875 

3000 

2319 

3055 

2895 

Pressure 
(mm  Hg) 

5.21 

5.21 

10.95 

10.95 

18.2 

Density 

<P/Po> 

9.84X10"4 

5.38X10"4 

1.68X10"3 

1.12X10-3 

2.19X10'3 

VDa 

98.5 

90 

98 

93 

96 
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TABLE  4.3.-  CONTINUED, 
(b)  Concluded. 


N.  Case 

No. 

36T 

37L 

37T 

38L 

38T 

39L 

Parameters^''*. 

Altitude 

(ft) 

200,000 

200,000 

200,000 

200,000 

200,000 

200,000 

Velocity 

(ft/sec) 

25,000 

10,000 

10,000 

10,000 

10,000 

10,000 

Body  half  angle 
(deg) 

20 

10 

10 

15 

15 

20 

Vq  (cir/sec) 

7.63x10s 

3.05x10s 

3.05x10s 

3.05x10s 

3.05x10s 

3.05x10s 

D  (cm2/ sec) 

4.4x10s 

3.2X104 

2.5X104 

1.7X104 

1.25X104 

l.oxio4 

b  (cm-1 ) 

3.42 

0.543 

1.785 

0.770 

3.25 

0.922 

Bq  (cm-1) 

173 

9.54 

12.2 

17.9 

24.4 

30.5 

Z 

1.16 

1.00 

1.00 

1.00 

1.00 

1.00 

Temperature 

(°K) 

3220 

1872 

1373 

1778 

1340 

1655 

Pressure 
(mm  Hg) 

18.2 

1.13 

1.13 

2.11 

2.11 

3.3 

Density 

(P/Po* 

1.73X10-3 

2.17X10-4 

2.94X10-4 

4.28X10-4 

5. 68X1 0"4 

7.34X10-4 

VDa 

93 

94 

95 

96 

97 

98 

4-23 


TABLE  4.3.-  CONCLUDED. 

(c)  Special  conical  geometry  cases  for  investigation 
of  effects  of  velocity  term. 


s.  Case 

Nv  No. 

13LA 

13LB 

13LC 

13LD 

13LE 

Parameters^"*. 

Altitude 

(ft) 

200,000 

200,000 

200,000 

200,000 

200,000 

Velocity 

(ft/sec) 

35,000 

35 ,000 

35,000 

35,000 

35,000 

Body  half  angle 
(deg) 

10 

10 

10 

10 

10 

VQ  (cm/ sec) 

1.068x10 

1.068x10 

1.068x10 

1.068X10 

1.068X10 

D  (cm2/sec) 

5.5X10® 

5.5X10® 

5.5x10® 

5.5x10® 

5.5x10® 

b  (cm-1 ) 

1.86 

0.54 

5.40 

54.0 

0 

Bq  (cm-1) 

194.5 

194.5 

194.5 

194.5 

194.5 

Z 

1.05 

1.05 

1.05 

1.05 

1.05 

Temperature 

(°K) 

2090 

2090 

2090 

2090 

2090 

Pressure 
(mm  Hg) 

8.65 

8.65 

8.65 

8.65 

8.65 

Density 

(P/PQ) 

1.45X10-3 

1.45X10-3 

1.45X10-3 

1.45X10-3 

1.45X10-3 

VD» 

99.5 

99.5 

99.5 

99.5 

99.5 

Velocity,  ft/sec 


Density,  lb/ft 


Enthalpy,  Btiylb 
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Figure  4.4.-  Boundary- layer  enthalpy  profile  at  5.76  feet 

from  conical  apex. 
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Figure  4.5.-  Profile  of  weak  field  collision  frequency  through 
the  boundary  layer  at  5.76  feet  from  conical  apex. 
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Distance  normal  to  surface,  ft 


Figure  4.6.-  Profile  of  strong  field  collision  frequency  through 
the  boundary  layer  at  5.76  feet  from  conical  apex. 
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5.  VARIATIONAL  METHOD  FOR  COMPUTATION  OF  THE  C-W  BREAKDOWN 

FIELD  STRENGTH 

In  this  section  a  variational  technique  will  be  fully  devel¬ 
oped  in  terms  of  orthogonal  polynomials  for  the  computation  of 
breakdown  field  strengths.  This  variational  principle  is  a  spe¬ 
cial  case  of  one  listed  in  Reference  26,  page  1109,  (see  also 
Appendix  B) . 

Suppose  Equation  (2.13)  is  multiplied  through  by  $  and 
integrated  over  all  space.  Upon  solving  for  the  eigenvalue  one 
obtains 


k2 


(5.1) 


where 


If  this  expression  vanishes  to  first  order  with  variations  in 
and  ^  (independently) ,  and  if  and  f  both  vanish  on  the 
boundaries  of  the  region,  then  if/  satisfies  Equation  (2.13)  and 
If/  satisfies  (see  Appendix  B) 

(V2  +  -  -  k2G$  (5.2) 

Here  $  is  called  the  adjoint  of  f  and  Equation  (5.2)  is  the 
adjoint  equation.  This  variational  principle,  6(k2)  ■  0,  is 
particularly  good  for  the  determination  of  eigenvalues  since  the 
error  in  k2  is  second  order  in  its  dependence  on  the  errors  in 
f  and  f  because  we  are  demanding  that  the  first-order  variation 
vanish. 

The  numerical  determination  of  k2  for  boundary-layer  flow 
will  now  be  formulated  for  several  cases  of  increasing  complexity. 
First,  the  special  case  of  constant  D  is  considered.  This 
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assumes  that  D  does  not  vary  significantly  in  the  region  of  appre¬ 
ciable  induced  electron  production.  This  approximation  should 
approximate  reality  fairly  well  since  Equations  (3.33)  and  (3.40) 

indicate  D  =  mc2/3mv  and  analysis  indicates  that  the  electron 

en 

RVDF  and  thus  (1/2 )mc2  varies  quite  slowly  with  electric  field 
strength  when  the  effective  field  strength  is  above  a  few  volts 
(Sections  3.3  and  3.4).  In  regions  of  lower  field  strength,  the 
RVDF  varies  considerably.  However,  there  is  no  increase  in  ioniza¬ 
tion  in  such  regions  and  diffusion  has  little  effect,  except  on 
the  downstream  side  of  the  applied  field.  Downstream  of  the  applied 
field  the  electrons  quickly  (in  about  three  mean  free  paths)  achieve 
a  Maxwellian  RVDF  among  themselves  but  remain  fairly  well  insulated 
from  the  heavy  particles.  Thus,  downstream  of  the  field,  the  aver¬ 
age  energy  will  change  considerably  (probably  to  about  one-half  its 
value  in  the  field) .  This  will  cause  some  error  in  the  downstream 
diffusion  rate  and  thus  in  the  eigenvalue  and  the  breakdown  field 
strength.  Likewise  there  is  a  corresponding  error  in  the  cross- 
flow  diffusion  rate  perpendicular  to  the  surface.  This  model  is 
formulated  below  in  two  and  three  dimensions. 

Following  the  constant-D  formulations,  the  variable-D  problem 
is  formulated.  This  latter  development  was  not  used  in  the  computa¬ 
tions,  however,  because  it  would  have  necessitated  the  use  of  either 
a  larger  computing  machine  or  a  significantly  slower  program.  The 
accuracy  of  the  input  data  for  the  computation  does  not  appear  to 
warrant  such  meticulosity . 

5.1  Formulation  for  Diffusion  Coefficient  Constant  Throughout 
the  Field 

When  D  is  constant.  Equations  (2.13)  and  (5.2)  become  iden¬ 
tical  except  for  the  sign  of  the  "b  •  V  term.  Since  is  in  the 
z  direction,  a  change  of  variable  z  -»  -z  transforms  one  equation 
into  the  other.  Thus  ^(-z)  -  ^(z),  so  that  f  need  not  be  varied 
as  an  independent  function.  It  will  be  seen,  subsequently,  that 
this  simplification  cuts  the  number  of  equations  to  be  considered 
exactly  in  half. 
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Let  us  assume  that  the  dependence  of  £  on  field  strength 
may  be  written 


(5.3) 


where  p*  ■  pZ(RT)CQl<j  “  (p/pQ) 2 (760)mm  and  p/pQ  is  the  density 
ratio,  pQ  being  the  density  under  STP  and  Z  the  compressibility. 
For  cold  air  this  polynomial  is  given  by  Equation  (3.49).  Then, 
with  Equations  (2.11)  and  (2.12),  G  becomes 


d.  d.+2 


(5.4a) 


and 


’cold 


- 37  s  I*iF 

^ai(El/p*1) 

5  F*  [F(f^)  -  30F]  l-rjof  <*•«» 


where 


di 

V&rrJ  aip  ,  _ 

ri  — d—  and  p-(iH- 

L  V 


760Z  (p  '/ p  ) 


mm-cm 

volt 


(5.5) 


and  for  truly  cold  air  Pcold  ■  p*.  Also,  from  Equation  (2.10) 


and  for  cold  air 


cold 


/  i  N2  ^  -  (d  .  +  2) 

■  (F* )  Iv 


-  (p*1)*  (0.243)  ^  ~330P^  xlO"* 


(5.6a) 


(5.6b) 


Thus  for  a  given  p* ' ,  P  fixes  k2  and  k2  fixes  P.  In  the 
following  formulation  we  have  initially  neglected  the  variation 
of  p*  and  have  taken  p*  ■  p*'  ■  value  of  p*  at  one-half  power 
point  (see  Section  5.3). 
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For  B  we  use  Equation  (4.5).  A  definite  electric  field  will 
now  be  chosen. 

5.1.1  Exponential  -  Gaussian  field  in  two  dimensions 

For  the  two-dimensional  problem,  it  is  assumed  that  all  quanti¬ 
ties  are  constant  in  the  x  direction  (flow  over  an  infinite  sur¬ 
face  and  an  infinite  antenna) .  This  problem  approximates  the  solu¬ 
tion  for  the  case  where  the  antenna  illuminates  the  gas  evenly  for 
a  very  long  distance  across  the  flow  as  compared  to  the  downstream 
illumination  width. 

Let  us  take  a  field  such  that  (see  Fig.  2.1) 

F  ■  E/E'  ■  exp  e^y  -  e2Z2^  (5.7) 

The  y-dependence  resembles  well  the  induction  field  measured  just 
outside  slot  antennas  (see  Section  6.1).  The  z-dependence  seems 
reasonable  but  is  chosen  arbitrarily.  Any  field,  however,  may  be 
represented  by  multiplying  Equation  (5.7)  by  polynomials  (e.g., 
Laguerre  y-representation  and  Hermite  z-representation) .  Also,  as 
will  be  evident  in  the  analysis  below,  little  extra  complication 
arises  with  such  representation  over  that  involved  from  the  use  of 
Equation  (5.7)  alone. 

For  a  trial  function  we  take  an  orthogonal  polynomial  expansion 

for 

^  -  Y  ir  „  (5.8) 

/  .  xn  i  n  in )  n 

min 

where  , 


‘  -  §  y 

~  _  2  z2  -  £“ 

i>  ■ 

rm,n 

e  VL°’2(y) 

e  2  H^(z)7  2 

- 

—  _ 

Here  L^’2  is  the  Laguerre  polynomial  given  by7 


7 

L„2  is  the  L  (y,a,3)  of  Reference  27. 
m  m 
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,a’2  -  eQy  y“2  ^IL  m+2'N 

m  dym  \  / 


(5.10a) 


and  is  the  Hermite  polynomial  given  by 


HI  -  (-l)n  .'l'2" 


n 


dz 


n 


(5.10b) 


The  A  are  the  parameters  to  be  varied,  and  y  and  a  are 
m,  n 

parameters  which  may  be  arbitrarily  selected  for  the  computation. 
However,  very  large  or  very  small  values  of  a  and  y  compared 
to  and  ,  respectively,  will  necessitate  the  use  of  more 

terms  in  Equation  (5.8)  for  good  accuracy. 

If  Equations  (5.8),  (4.5),  and  (5.4a)  or  (5.4b)  are  substi¬ 
tuted  into  Equation  (5.1)  with  ^  ■  if/(- z)  ,  one  obtains  (Appen¬ 
dixes  C,  D,  and  E) 


-  k2 


I  I  I-1’"' 

m,n  m' ,n' _ 


I  I 

m,n  m' ,n' 


(-I)”'  *0".  A  ,  i  A 

m  ,n  in  >n  ro^n 


(5.11) 


where 


n' 


+  nfn-D-y2! 


Y,o 

n '  ,n 


+  I 


Y,o 
n'  ,n 


(nrt-3) 


Jm'°m  "  a  2) 

m  ,m 


ja;°  +  al  ja,o 
m  ,m  4  m'  ,m 


(5.12) 


is  the  H  (z,a)  defined  in  Reference  27. 
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Nn'»n  ■  b  y 
m*  ,n'  Bo  ’ 


-C**  -  j°»b  ^  C-  i  lY>° 

m',m  um',my  ^  2  V, 


°  +  nv  i*»° 

n+i  nY  An'  ,n 


Mm,n  -  v 
m '  ,  n '  Y 


I  V 


a >  (dj+2 ) e 

J  ,  1  1 

m  ,m 


7, (d.+2) e 

1  c 


In 1  ,n 


or,  for  cold  air, 


(*"-). 


,-w  x 


cold 


1  -  30P 


Here 


and 


*,*e\  /Y.3€\ 
m' ,m J  V  ‘ ln' ,n  J 


~  30pl  Jm* 


a, 26^/7, 2€X 
,m  i^n’.n  J 


00 

i™  .  f  H“  H=  dz  .  jY.j 

m,n  /  m  n  n,m 


c  w 

j01’^  -  f  e-(a+n)y  yc  l«»2  l“»2  dy  - 
m,n  J  1  m  n  *  n,m 


It  is  shown  in  Appendix  D  that 


t  Y»n 
n,n 


2 (2n-l) 


n-i 


4 (n-1)2 


yIIx-ji) 

Tfrj 


xY.n 
n-2 ,  n-s 


TY»n  .  /'njtlv.'W  j7,ti  1  _ ( H±jA  t7 » q 

m,m+-2V  \mf  1  /  t}  m+2v-i  ,nH-i  2  (nfl-1)  \  7T)  J  nH-2V,mt2 

and  I  ■  0.  The  lowest  members  are 

m,m+2V+i 


-0 

(5.13) 

(5.14a) 

(5.14b) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 
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and 


I 


Y.T) 

0,0 


rY.n  .  2y2  /  Tr 

1,1  7  +  n  v  y  +  n 


Similarly  in  Appendix  E  it  is  shown  that 


<“+1»  «  ■  [(n-»l-c)«(2«Dn] 


+  m(mf2)ot 

m-i j  n— i 


(5.19) 


The  lowest  members  are 


c 

ja.T 


cl 


(a  +  t|) 


-  , 

■1  1,0 


cl 


(a  +  tj) 


C+  2 


|^(2  -  c)a  +  3t]J 


These  reduction  relations  allow  all  the  I's  and  J's  to  be 
generated  from  the  lowest  members.  The  normalization  integrals 
for  Hermite  and  Laguerre  polynomials  are  special  cases  of  the  I's 
and  J's.  They  are 


I 


7,o 

m,n 


Y^"n:(27)n 


6n 

m 


^a,o  _  n! (n+2) I  &n 
m,n  q3  m 


where  6n  is  a  Kronecker  delta, 
m 

If  one  now  evaluates  dks/dA  from  Equation  (5.11)  and  sets 

m,n  ^ 

them  equal  to  zero,  as  required  by  the  variational  principle,  one 
finds9 


Here  the  symmetry  of 


of 


Nm’n 
m ' ,  n ' 


<C?n' 


and  M^|nn,  and  the  antisymmetry 


in  the  primed  and  unprimed  indices  has  been  used, 


for  example ,  Q 


m ,  n 
m '  ,  n ' 


-  Q^'  >n' 

Tn,n 


but 


^n,n 
m ' ,  n ' 


■  _  ' »n' 

Mm,n 
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Setting  the  determinant  of  coefficients  of  this  set  of  equations 

(for  the  A  )  equal  to  zero  determines  k2 .  Alternately,  one 
m ,  n 

could  specify  one  of  the  A  and  use  all  but  one  of  these  equa- 

m ,  n 

tions  along  with  Equation  (5.11)  to  determine  k2 .  These  approaches 
are  completely  equivalent.  We  want  the  lowest  value  of  k2 ,  that 
is,  we  want  the  lowest  value  of  k2  that  permits  the  vanishing  of 
the  determinant. 

5.1.2  Exponential  -  Gaussian  field  in  three  dimensions 

The  extension  of  Section  5.1.1  to  three  dimensions  is  straight¬ 
forward.  For  the  electric  field  we  now  take 


-e  y  — €  z2-e  x2 
F  -  E/E'  -  e  1  e  2  3 


(5.21) 


Again,  as  discussed  earlier  after  Equation  (5.7),  multiplication 
by  various  polynomials  allows  adequate  representation  of  any  field 
without  introducing  significant  extra  work. 

In  place  of  Equation  (5.8)  we  have 


l 


A 


l ,  m ,  n 


*1 


,m,  n 


(5.22) 


with 


a 

-  7  y 


l  ,m,n 


l£’2(y)e 


(ex2+-yz2) 


H 


(x)H^(z)p 


1 

2 


y 


n, 

2 


(5.23) 
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and  for  cold  air 


M 


l 

l 


,m,n 

I  m  I  —I 

,m  ,n 


-  |-(n+n')  -  i(l+l') 

y  3 


2 

a,3€ 

J  ,  _ 
m  ,m 


Y»3€_ 

V,n 


-  30P  I 


3  |2£, 
IM* 


2 

CX|2€1 
Jm'  ,m 


(5.27b) 


5.2  Formulation  for  Nonconstant  Diffusion  Coefficient 

In  Sections  3.3  and  3.4  it  was  pointed  out  that  the  diffusion 
coefficient  is  a  function  largely  of  the  local  composition  in  the 
regions  of  interest  and  of  the  heavy  particle  temperature.  In  addi¬ 
tion,  the  diffusion  coefficient  depends  on  the  electric  field 
strength  also  in  the  weak  field  regions,  that  is,  field  strengths 
less  than  about  10  volts.  Although,  as  pointed  out  earlier,  one 
can  get  fairly  good  results  with  the  constant-D  formulation,  a 
truly  accurate  treatment  for  the  case  of  strong  mixture  gradients 
and  large  changes  in  the  electric  field,  both  of  which  are  absent 
in  the  usual  closed  cavity  breakdown  studies ,  should  include 
variable-D  effects. 

In  general,  it  will  always  be  possible  to  obtain  reasonable 
fits  to  data  in  the  form 


T 


I 

i>  j  >k 


-V 


i»  j  »k 


E 


c 


i>  j 


(5.28) 


(5.29) 


When  D  is  not  constant,  one  cannot  obtain  ip  from  ip  by 
merely  taking  ip{ z  -►  -z).  The  trial  function  for  f  must  be  varied 
independently.  The  trial  function  for  ip  may  again  be  taken  as 
Equation  (5.8)  with  (5.9)  for  the  two-dimensional  case  and  (5.22) 
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with  (5.23)  for  the  three-dimensional  case.  Since  the  three- 
dimensional  case  is  a  simple  generalization  of  the  two-dimensional 
case,  only  the  development  of  the  two-dimensional  case  will  be 
fully  demonstrated. 

5.2.1  Exponential  -  Gaussian  field  in  two  dimensions, 
nonconstant  D 

For  the  $  trial  function  we  may  take 

?  -  V  B  J  (5.30) 

/  .  m,n  in f n 

m»n 


Here  again  ^  is  given  by  Equation  (5.9). 
m  i  n 

For  k2  we  find 

E  E  (fen'  -  *0?n')  V ,»•**,» 

_  k2  .  P.i.n.jn.'.x.n,: _ 

E  E  V.n'Vn 

m,n  m' ,n' 


(5.31) 


Here  Q™ J n„  •  is  still  given  by  Equation  (5.12);  however,  N1"?"  , 
m  ,  n  m  ,  n 

and  M™|nn,  are  now  changed.  For  any  specific  case  they  may  be 

readily  developed  by  multiplying  Equations  (5.28)  and  (5.29)  by 

V*  i  and  .  integrating  and  making  use  of  the  reduction 

m,n  m,n 

formula  (C.2a)  (see  Appendix  C)  and  the  integrals  (5.15)  and  (5.16). 

The  set  of  equations  dk2/dA  ■  dk2/dB  •  0  become 

m ,  n  m ,  n 


X'  fen'  - 

+  k2 

o 

■ 

c 

at 

«E 

m'  ,n' 

/  ( Qm’n  +  Nra,n 

l—j  \  m'  ,n'  m'  ,n' 

*  k* 

)  V.n'  -  0 

m' ,n '  N 

S 

where  we  have  made  use  of  the  symmetry  of  Q  and  M  and  anti¬ 
symmetry  of  N  with  respect  to  the  primed  and  unprimed  indices. 
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5.3  Discussion  of  Approximation  that  Particle  Density  or  p* 
is  Constant 

At  the  end  of  the  introduction  to  Section  5.1  the  approxima¬ 
tion  p*  ■  p*‘  was  made  and  the  value  selected  for  these  quanti¬ 
ties  was  the  value  of  p*  at  a  distance  above  the  aperture  center 
corresponding  to  one-half  the  field  strength  at  the  aperture  (see 
Section  7.1).  Since  p*  ■  760Z(p/pQ),  this  is  equivalent  to 
assuming  that  the  particle  density  is  constant  in  the  boundary 
layer.  The  error  involved  should  not  be  too  important  in  view 
of  the  uncertainty  in  the  values  already  used  for  the  ionization 
frequency. 

The  consideration  of  the  variable-p*  problem  is  actually  no 
more  difficult  than  the  constant-p*  problem.  This  may  be  done 
by  using  a  polynomial  fit  to  the  y-dependence  of  p*.  The  addi¬ 
tional  integrals  occurring  in  the  variational  principle  give 

extra  terms  in  Equation  (5.14b)  which  may  be  written  down  immeai- 

c 

ately  in  terms  of  the  and  's.  This  unnecessary 

approximation  was  actually  the  result  of  an  oversight.  Since 
almost  all  experimental  breakdown  studies  have  been  conducted  at 
room  temperature,  the  macroscopic  coefficients  derived  from  these 
experiments  have  usually  been  represented  as  functions  of  E/p, 
where  p  was  the  actual  gas  pressure,  although  the  truly  meaning¬ 
ful  physical  variable  should  be  E  divided  by  particle  density. 

In  our  initial  calculation,  therefore,  we  had  used  p  instead  of 
p*  in  Equations  (5.4)  to  (5.6).  Also,  p  is  constant  through 
the  boundary  layer  so  that  p  -  p' .  Thus,  effectively,  the  compu¬ 
tations  were  conducted  with  the  net  ionization  frequency  corre¬ 
sponding  to  a  constant  density  of  about  a  factor  of  three  to  an 
order  of  magnitude  too  high  but  with  the  diffusion  constant  at 
the  proper  density.  These  results  are  presented  in  Section  7.2.1, 
although  incorrect,  because  they  depict  qualitatively  all  the 
phenomena  to  be  expected  with  breakdown  in  a  boundary  layer .  When 
the  oversight  was  noted,  some  results  were  quickly  recomputed  at 
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the  more  reasonable,  but  still  constant,  particle  densities  and 
the  results  presented  in  Section  7.2.2. 

There  is,  of  course,  no  constant-p*  approximation  necessary 
in  the  zero-velocity  (no  flow)  breakdown  computations  of  Section  6. 
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6.  ZERO  VELOCITY  BREAKDOWN  COMPUTATIONS  AND  COMPARISON  WITH 
EXPERIMENT 

In  this  section,  computations  for  zero  velocity  will  be  com¬ 
pared  with  experiments  conducted  with  antenna  configurations 
similar  or  equivalent  to  a  type  of  slot  antenna  commonly  used  for 
reentry  vehicles  (Ref.  4). 

6.1  The  Shape  of  the  Electric  Field 

The  breakdown  studies  of  Reference  4  includes  probe  measure¬ 
ments  of  the  electric  field  strength  of  slot  antennas  at  various 
distances  away  from  the  antenna  surface  along  the  line  perpendicular 
to  the  aperture  surface  of  its  geometric  center.  These  data  are 

plotted  in  Figure  6.1  as  E/E  versus  y/d,  where  y  is  distance 

s 

above  the  aperture  surface,  d  is  the  small  dimension  of  the 

aperture,  and  E  is  the  field  strength  at  the  aperture  surface, 
s 

Notice  that  all  the  data  cluster  quite  well  about  the  representative 
curve  drawn  except  for  one  of  the  flared  horns. 

Although,  as  pointed  out  in  Section  5.1.1,  it  is  possible  to 
match  this  representative  curve  with  any  desired  accuracy  in  the 
breakdown  computations,  it  was  found  sufficient  to  merely  match 
the  initial  slope.  Thus,  in  Equation  (5.7)  is  taken  to  be 

€  -  0.686/d  (6.1) 

For  the  electric  field  variation  across  the  gap,  there  are 
no  data.  A  theoretical  computation  of  such  a  near  field  is  very 

difficult.  We  have  taken  e  in  Equation  (5.7)  such  that  the 

2 

power  falls  to  half  at  z  m  +_  d/2,  that  is, 


or 

.  1.388 
2  d2 


(6.2) 
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This  type  of  field  variation  will  reasonably  represent  the  "size" 
of  the  field  and  seems  to  yield  good  agreement  between  computed 
and  experimental  breakdown  field  strengths. 

The  x-dependence  is  not  important  here  because  the  aperture 
is  taken  to  be  very  long  in  the  x-direction  so  that  a  two-dimensional 
treatment  applies  (see  Section  5.1). 

The  relationship  between  breakdown  field  strength  and  the  radi¬ 
ated  power  level  at  breakdown  may  be  obtained  using  a  result  from 
Reference  4,  page  81.  In  this  reference  it  is  shown  that  for  a 
long  narrow  slot,  less  than  1/2  wavelength  long  the  power  is  related 
to  the  slot  field  strength  by 

P  -  E  2d2  —  (6.3) 

S  „  2 


where  d  is  the  slot  width,  z  is  the  complex  slot  impendance, 

s 

and  R  is  the  real  part  of  z  .  In  turn,  z  is  related  to 
s  s  s 

zd,  the  impedance  of  a  dipole  of  diameter  d/2.  This  relation  is 


z 

s 


70,952 

Zd 


(6.4) 


The  dipole  impedance  may  be  determined  from  Reference  32  (pp.  154-179) . 

6.2  Convergence  with  4x4  Term  Representation 

The  coefficient  functions  for  Equation  (5.20)  were  generated 
as  indicated  in  Section  5.1.  The  lowest  eigenvalues  of  Equation 
(5.20)  were  determined  by  solving  for  the  k2  that  allows  the 
determinant  of  Equation  (5.20)  to  vanish.  This  was  done  by  itera¬ 
tion  on  an  IBM  1620.  The  breakdown  voltage  is  then  given  by  Equa¬ 
tions  (5.5)  and  (5.6b).  Our  computations  have  utilized  four 
Hermite  polynomials  and  four  Laguerre  polynomials.  This  gives  a 
total  of  sixteen  coefficients  to  vary  or  a  16x16  determinant. 

Theory  tells  us  (see  Ref.  26  regarding  variational  techniques 
with  positive  definite  operators)  that  the  better  the  accuracy, 
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the  lower  the  computed  eigenvalue;  that  is,  the  eigenvalue  is 
approached  from  above.  Zn  the  present  case,  this  implies  that  any 
alteration  in  the  trial  functions  that  produces  a  lower  breakdown 
voltage  is  a  desirable  alteration  and  gives  a  more  accurate  eigen¬ 
value  . 

To  examine  the  accuracy  of  the  4x4  term  representation,  the 
breakdown  field  was  computed  with  2x2,  3x3,  and  4x4  terms  with  var¬ 
ious  values  of  a  and  y  and  the  results  plotted  in  Figure  6.2. 

It  is  seen  that  with  a  “  2ei  and  y  -  1/2 the  change  of  the 
breakdown  field  is  not  large  with  addition  of  terms  and  that  the 
rate  of  change  at  n  ■  4x4  is  quite  small.  This  indicates  that 
our  results  are  at  most  2  volts,  or  less  than  5  percent,  above  the 
result  which  would  be  obtained  for  this  model  using  many  more 
terms  in  the  expansion.  All  our  computations  used  a  ■  2e^  and 

y  -  1/2 e  . 

'  -  2 

It  should  be  noted  that  2x2  representation  is  so  crude  that 
it  does  not  even  allow  variation  of  the  z-dependence  curvature  in 
the  electron  density  function  at  z  ■  0. 

6.3  Comparison  of  Breakdown  Field  Strength  with  Experiment 

The  only  controlled  experiments  on  antenna  breakdown  in  air 
have  been  performed  in  ambient  temperature  air  at  rest.  The  results 
(Ref.  4)  of  one  series  of  these  experimental  measurements  for  a 
17-inch  long  slot  antenna  of  width  d  operating  at  380  megacycles 
is  given  in  Figure  6.3  (solid  line).  For  comparison  purposes,  the 
results  of  a  two-dimensional  constant-D  model  computation  of  break¬ 
down  are  given  as  a  dotted  line.  For  these  calculations,  the  values 
y  ■  O.Se^  and  a  ■  2.0^  were  used  with  slot  widths  of  0.1,  1.27, 
and  2.54  centimeters. 

The  computed  points  agree  well  with  experiment  in  the  region 
in  which  the  input  data  for  £  are  valid,  that  is,  for  E/p  110 
volts-cm  1-mm-1Hg.  Marked  deviation  from  the  experimental  curve 


k 


6-4 


occurs,  as  should  be  expected,  for  E/p  >  200  volts -cm-1 -mm"1 , 
since  the  extrapolation  of  measured  macroscopic  coefficients  into 
this  regime  is  expected  to  yield  poor  values  for  The  com¬ 

parison  supports  the  validity  of  the  model  for  the  simple  geometry 
investigated  here. 


Figure  6.1.-  Variation  of  field  strength  with  distance  for 


Value  of  eigenvalue,  volts/ cm- 


E  /p ,  volts/cm- 


pd ,  nun  Hg-cms 

Figure  6.3.-  Comparison  of  breakdown  computation  with  experl 
ment  for  a  narrow  slot  and  zero  velocity  flow. 
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7.  COMPUTATIONS  WITH  FLOW  RESEMBLING  BOUNDARY  LAYERS  ENCOUNTERED 
IN  HIGH-SPEED  FLIGHT 


7.1  Review  of  Input  Assumptions 

The  inviscid  and  viscous  layer  properties  for  cones  and 
wedges  at  various  altitudes  and  reentry  velocities  were  computed 
as  described  in  Section  4.  The  boundary-layer  velocity  profile 
divided  by  D  was  represented  as  discussed  in  Section  4.3.  In 
this  representation  of  V/D,  it  was  assumed  that  the  diffusion 
was  fully  ambipolar  and  the  diffusion  coefficient,  computed  for 
hot  air  as  described  in  Section  3.4.1,  was  evaluated  for  gas  tem¬ 
peratures  and  densities  at  a  distance  from  the  antenna  surface 
corresponding  to  a  point  where  the  electric  field  intensity  is 
one  half  of  its  value  in  the  aperture.  The  reason  for  this  choice 
is  that  the  ambipolar  theory  requires  that  D  be  evaluated  at 

cl 

the  point  where  breakdown  first  tends  to  occur.  This  would,  for 
example,  be  at  the  center  of  a  cavity  for  closed  cavity  breakdown. 

.  .  „  .  j-1/2 

It  would  be  about  a  diffusion  length 


,  A  -  ^  d2  (nD)/ dy2J  / (nD)J 


from  the  surface  in  the  present  case  and  would  require  iteration 
for  A  to  evaluate  properly.10  The  error  will  not  be  large  by 
the  choice  indicated.  The  more  accurate  iteration  treatment  is 
not  warranted  because  of  the  poor  state  of  knowledge  of  the  data 
used  for  the  other  inputs. 

The  value  of  the  ionization  coefficient,  £ ,  is  taken  as 
indicated  by  the  choice  of  v  in  Section  3.4.2.  It  is  given  by 
Equation  (3.54) 


^hot 


^cold(p) 


10That  is,  select  A,  evaluate  D_  at  this  point, 

a 


^  £d2 (nD)/dy2 


-i/s 


» 


compare,  correct  A, 


solve  for 
and  repeat. 
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where  £col<3  is  given  by  Equation  (3.49)  and  P  «  p*/E.  This 

value  of  C.  .  is  taken  such  that  D  was  evaluated  a  distance 
hot  a 

from  the  aperture  as  discussed  above  for  V/D  .  and  p*' ,  the 
effective  pressure  as  defined  on  page  3-26,  taken  at  the  center 
of  the  antenna  aperture,  was  evaluated  as  discussed  in  Section  5.3. 
Initially,  as  noted  in  Section  5.3,  due  to  an  oversight,  many 
computations  were  made  with  the  ionization  frequency  evaluated  at 
an  improperly  high  density  corresponding  to  p  rather  than  p*: 

This  is  equivalent  to  having  a  net  ionization  frequency  of  about 
an  order  of  magnitude  or  more  higher  than  the  actual  values. 

Although  incorrect,  these  results  do  depict  qualitatively  the  type 
of  breakdown  behavior  expected  in  a  boundary  layer  and  are  there¬ 
fore  discussed  in  Section  7.2.1.  All  such  results  are  clearly 
labeled  illustrative  only.  With  the  oversight  corrected  some  of 
the  computations  were  repeated  and  these  results  are  discussed 
in  Section  7.2.2. 

The  electric  field  adopted  approximates  that  for  a  narrow 
slot  antenna,  as  described  in  Section  6.1,  and  corresponds  to  a 
gap  width  of  1/2  inch.  The  arbitrary  parameters  a  and  7  were 
taken  as  in  Section  6.2  for  good  convergence  with  4x4  repre¬ 
sentation  when  there  is  no  flow  (a  ■  2^,  7  ■  (1/2)6^).  This 

model  compares  well  with  experiment  for  cold  air  without  flow  in 
the  same  geometry  (Section  6.3).  However,  it  was  found  that  with 
flow  the  4x4  term  representation  used  proved  inadequate  for 
many  of  the  computations  of  Section  7.2.2  with  the  proper  density 
or  p*  used,  and  further  results  would  require  the  use  of  more 
terms.  The  method  of  solution  is  the  same  as  in  Section  6, 
namely,  the  determinant  of  the  coefficients  of  Equation  (5.20)  is 
set  equal  to  zero  and  solved  by  iteration  for  k2  or  P,  the 
relationship  between  k2  and  P  being  given  by  Equation  (5.6b). 

Since  the  gap  antenna  was  assumed  very  long,  the  two-dimensional 
formulation  of  Section  5.1.1  was  used  for  both  cones  and  wedges. 

The  development  of  the  variational  principle  used  could  be  reformu¬ 
lated  in  cylindrical  coordinates  for  antennas  mounted  in  cones. 
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In  fact,  there  would  be  fewer  terms  in  such  a  formulation  because 
the  Laguerre  polynomials  behave  more  simply  with  the  cylindrical 
coordinate  differential  operator.  However,  the  antenna  is  assumed 
to  be  placed  sufficiently  far  downstream  in  this  model  that  the 
conical  problem  may  quite  accurately  be  approximated  by  the  "two- 
dimensional11  Cartesian  problem  and,  therefore,  a  separate  program 
was  unwarranted. 

Computations  were  attempted  for  several  flight  velocities 
between  10,000  and  45,000  feet  per  second  at  altitudes  of  200,000 
and  250,000  feet  for  both  cones  and  wedges  of  10°,  15°,  and  20° 
half  angles.  The  gap  antenna  was  assumed  placed  at  a  distance 
sufficiently  downstream  that  a  2-foot  body  total  width  or  diameter 
existed  at  that  point.  This  was 

*■<*>  -  /(5?JJIS2)  ‘’-1' 

At  200,000  feet  these  conditions  made  transition  to  turbulent 
flow  probable,  so  that  computations  for  this  altitude  were  attempted 
for  both  a  laminar  and  a  turbulent  boundary  layer.  However,  even 
with  use  of  a  proper  ionization  frequency,  the  results  of  the  calcu¬ 
lations  for  the  turbulent  cases  cannot  be  considered  adequate  since 
the  effective  electron  diffusion  rates  are  probably  significantly 
different  in  a  turbulent  region  (see,  e.g.,  Ref.  35),  whereas  they 
are  assumed  unaltered  by  turbulence  in  this  computation. 

7.2  Discussion  of  Computed  Breakdown  Field  Strengths 

7.2.1  Illustrative  computations  using  improper  ionization  frequency 

The  input  data  for  the  illustrative  computations  conducted 
with  the  ionization  frequency  evaluated  at  improper  particle  densi¬ 
ties  (discussed  in  Sections  7.1  and  5.3)  are  listed  in  Table  4.3, 
Cases  2  through  39,  and  the  results  of  these  computations  are 
depicted  in  Figures  7.1  through  7.6. 
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In  Figures  7.1  through  7.4  the  breakdown  field  strength11  in 
the  antenna  aperture  is  plotted  versus  the  vehicle  velocity.  It 
is  seen  that,  in  all  cases,  the  breakdown  field  strength  climbs 
significantly  with  vehicle  velocity  and  thus  this  appears  to  sug¬ 
gest  that,  other  things  being  equal,  one  may  transmit  at  a  sig¬ 
nificantly  higher  power  level  at  higher  velocities.  When  plotted 
as  Eg/p,  as  in  Figures  7.5  and  7.6,  the  curves  are  relatively 
flat  and  even  slightly  decrease  with  increasing  vehicle  velocity. 

This  implies  that  the  increase  in  computed  breakdown  field  strength 
with  velocity  is  largely  a  pressure  (or  actually  density  -  see 
Section  5.3)  effect,  that  is,  due  to  the  increased  shock-wave  com¬ 
pression  at  higher  velocities.  Notice  also  the  effect  of  density 
in  the  larger  values  of  Efi  at  200,000  feet  relative  to  those  at 
250,000  feet,  whereas  the  values  of  Eg/p  at  these  two  altitudes 
are  only  slightly  different. 

In  order  to  examine  the  effect  of  the  flow  velocity  profile 
itself  at  the  antenna  position,  Case  13L  (laminar)  was  rerun 
several  times  with  various  values  of  b,  all  other  parameters 
remaining  constant.  This  adjusts  the  velocity  profile  so  that  it 
climbs  to  the  free-stream  value  at  various  rates  with  distance 
away  from  the  wall .  Very  high  b  corresponds  to  a  boundary-layer 
velocity  profile  of  negligible  thickness  -  as  if  the  flow  were 
uniform  over  the  surface.  Smaller  values  correspond  to  thicker 
velocity  profiles,  whereas  b  ■  0  corresponds  to  no  flow  or, 
equivalently,  an  infinite  boundary-layer  thickness.  The  input 
data  for  these  calculations  are  listed  as  Cases  13L,  13LB,  13LC, 
and  13LD  of  Table  4.3,  and  the  results  are  plotted  in  Figure  7.7. 

It  is  seen  that  by  sweeping  out  electrons  the  local  flow  velocity 
does  play  a  significant  role  in  breakdown  at  the  actual  boundary-layer 

13"Note  that  the  quantity  determined  here  is  Ee,  the  effective 
field  strength  defined  by  Equation  (3.19a).  For  air 
E 

E  «  - - 7 —  .  For  p.  >  200  mm  Hg-cm,  E  -  E„  . 

e  r  -i  i/ 2  *  e  rms 

[l  -  (36/PX)2] 
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thickness  (but,  recall,  an  improperly  high  ionization  frequency) 
and  also  that  a  velocity  profile  with  b  ■  10e?  or  larger  will 
give  about  the  same  result  as  uniform  flow  over  the  surface. 
However,  the  density  effect,  due  to  shock-wave  compression,  is 
considerably  larger  than  the  effect  of  local  flow  velocity.  The 
velocity  effect  would  be  stronger  if  a  more  proper  ionization 
frequency  were  used  (Section  7.2.2). 

7.2.2  Computations  with  a  more  proper  ionization  frequency 

Computations  were  attempted  with  the  corrected  ionization 
frequency  for  the  same  vehicles,  vehicle  velocities,  and  altitudes 
as  were  used  in  Section  7.2.1.  However,  it  was  discovered  that 
only  a  very  few  of  the  eigenvalues  or  breakdown  voltages  could 
then  be  found  using  the  4x4  term  representation  (16  constants 
to  vary),  that  is,  with  4  terms  each  for  the  y-  and  z-expansions , 
respectively.  Analysis  shows  that  the  use  of  a  few  more  terms  in 
the  trial  function  would  alleviate  this  problem.  In  fact,  4x5 
term  representation  (5  for  the  z-dependence  at  each  y)  would 
probably  give  all  the  eigenvalues  but  4x7  is  desirable  for  good 
accuracy . 

Before  discussing  the  analysis  leading  to  these  conclusions 
and  presenting  the  few  computed  results  obtained,  let  us  examine 
why  this  difficulty  arose  here  but  not  in  Section  7.2.1,  where 
the  velocities  were  the  same  but  the  ionization  frequency  was 
evaluated  at  the  wrong  p*  or  density. 

If  the  quantity  p*,  the  modified  pressure  defined  in  con¬ 
nection  with  Equation  (5.3),  were  actually  utilized,  instead  of 
the  improper  p,  and  if  p*  is  taken  constant  over  all  space  (at 
(•gr)  ■'2)  in  Equations  (5.5)  and  (5.6),  the  effect  is  as  if  the 

net  ionization  frequency  were  a  smaller  function  of  P  than  used 
in  Section  7.2.1  by  the  factor  (p*/p)2*  This  factor  varies 
between  about  10_1  and  10“2 .  The  functional  dependence  of  the 
ionization  frequency  on  P  ■'*  p*/E*  and  F(y,z)  ■  E/E 1  is 
(the  product  of  Eqs.  (5.4b)  and  (5.6b)) 
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£  «  P2(P  -  30P)/P3  (7.2) 

When  the  coefficient  of  this  function  is  reduced  by  the  factor 
(p*/P)s»  the  quantity  P  must  become  smaller  at  breakdown. 

That  this  is  so  becomes  clear  if  one  notes  in  Equation  (7.2) 
that  when  P  is  smaller,  the  region  of  space  above  the  aperture 
over  which  v  is  positive,  and  thus  gives  net  electron  production 
rather  than  loss  by  attachment,  is  increased.  This  causes  the 
diffusion  to  be  spread  over  a  larger  volume  and  results  in  a 
decreased  diffusion  loss  rate.  This  decreased  diffusion  rate 
is  necessary  for  the  breakdown  conditions  since  the  ionization 
frequency  has  been  reduced  by  the  factor  (p*/p)2«  If  there  were 
no  flow,  both  the  overall  diffusion  term  and  net  production  term 
are  smaller,  but  balance  at  breakdown.  Thus  here  the  velocity 
convection  term  need  not  be  as  large  to  effect  the  breakdown 
relation  as  in  Section  7.2.1.  Since  the  velocity  term, 

(v/D)  •  V(nD),  tends  to  vary  inversely  as  the  length  of  the  field 
containing  region  whereas  the  diffusion  term  tends  to  vary  inversely 
as  the  square  of  such  a  length,  one  should  expect  phenomena  asso¬ 
ciated  with  local  flow  velocity  to  be  significantly  larger  due  to 
the  smaller  ionization  frequency.  This  is,  in  fact,  what  was  found. 

In  order  to  examine  numerically  the  behavior  of  these  eigen¬ 
values  with  increasing  free-stream  velocity,  the  input  data  of 
Table  4.3,  Case  5,  (with  the  more  proper  value  p*)  were  used 
except  with  varying  values  of  Bq.  With  the  proper  value  of 
Bq  ■  31.8,  no  eigenvalue  could  be  found  for  these  conditions  with 
4x4  term  representation.  With  the  use  of  an  even  more  crude 
2x2  term  representation,  the  determinant  of  the  coefficients 
of  Equation  (5.20)  was  plotted  versus  P  for  various  values  of 
BQ  in  Figure  7.8.  (Note  the  use  of  a  hybrid  scale  for  the  value 
of  the  determinant.)  The  first  crossing  of  the  "zero  axis"12 

18 It  should  be  noted  that  for  plotting  convenience  in  Figure  7.8 
a  value  of  the  determinant  of  +103  was  considered  small  enough, 
relatively,  to  be  considered  zero. 
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gives  the  eigenvalue  (in  terms  of  P  which  is  uniquely  related 
to  k2  and  E' )  for  each  case.  It  is  seen  that  the  eigenvalue 
disappears  between  BQ  ■  6  and  B  ■  7.  The  4  by  4  term  repre¬ 
sentation  was  then  used  with  BQ  ■  7  and  the  eigenvalue  readily 
computed  and  found  to  be  P  ■  3.33  or  E'  ■  22.3  volts/ cm.  This 
demonstrates  that  the  use  of  more  terms  does  allow  one  to  obtain 
eigenvalues  that  could  not  be  obtained  with  fewer  terms. 

The  "physical  reason"  for  the  above  phenomena  may  be  seen 
as  follows:  The  probable  shape  of  the  z-dependence  (at  a  con¬ 
stant  y)  of  the  eigenfunctions  are  sketched  schematically  in 
Figure  7.9  for  expansions  using  two-term,  three-term,  four-term, 
and  five-term  functions  corresponding  to  a  high  value  of  Bq. 

Since  H  ,  the  Hermite  polynomial,  is  a  polynomial  of  degree  n 
and  is  odd  if  n  is  odd  and  even  if  n  is  even,  these  eigenfunc¬ 
tions  behave  like  polynomials  of  first,  second,  third,  and  fourth 
order  multiplied  by  e~^z  ^2.  If  the  loops  in  Figure  7.9  contain¬ 
ing  a  region  of  negative  slope  to  the  left  of  the  actual  maximum 
(velocity  to  the  right)  occur  in  a  region  where  v  is  positive 
then  only  "breakdown"  can  occur  and  there  can  be  no  critical 
breakdown  condition.  This  is  because,  at  a  point  in  such  a  loop, 
the  diffusion  term  is  toward  the  point,  tending  to  increase  the 
electron  density,  and  in  the  forward  part  of  the  loop  (slope  down) 
the  velocity  term  is  also  tending  to  increase  the  electron  density. 
Thus  all  terms  in  the  conservation  equation  "produce"  electrons. 

The  phenomena  is  absent  at  lower  velocities  because  the  asymmetry 
demanded  of  the  eigenfunction  is  smaller  at  lower  velocities. 

Finally,  in  Figures  7.10  and  7.11  we  present  the  only  eigen¬ 
values  which  were  obtained  using  the  4  by  4  representation.  The 
field  strengths  here  may  be  somewhat  in  error  for  the  model  since 
high  accuracy  is  not  possible  here  using  only  a  4  by  4  term  repre¬ 
sentation.  The  short  lines  through  the  single  points  indicate  the 
estimated  behavior  of  the  breakdown  curve. 


Breakdown  field  strength,  volts/ cm 
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Figure  7.1.-  Aperture  breakdown  field  strength  versus  vehicle  velocity 
for  a  cone  at  200,000-foot  altitude  and  zero  angle  of  attack.  For 
purposes  of  illustration  only.  (Improper  ionization  frequency  used. 
See  Section  7.1.) 
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Vehicle  velocity,  thousands  of  feet  per  second 


Figure  7.2.-  Aperture  breakdown  field  strength  versus  vehicle  velocity 
for  a  wedge  at  200,000-foot  altitude  at  zero  angle  of  attack.  For 
purposes  of  illustration  only.  (Improper  ionization  frequency  used. 
See  Section  7.1.) 
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Figure  7.3.-  Aperture  breakdown  field  strength  versus  vehicle  velocity 
for  a  cone  at  250,000-foot  altitude  at  zero  angle  of  attack.  For 
purposes  of  illustration  only.  (Improper  ionization  frequency  used. 
See  Section  7.1.) 
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Figure  7.4.-  Aperture  breakdown  field  strength  versus  vehicle  velocity  for 
a  wedge  of  250 ,000-foot  altitude  at  zero  angle  of  attack.  For  purposes 
of  illustration  only.  (Improper  ionization  frequency  used.  See 
Section  7.1.) 
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Eg/p,  volts/ cm-i 


7-13 


80 

60 

40 

20 

r 

r 

o 

80 

60 

40 

20 

> 

0 

80 

60 

40 

20 

** 

0 


A. 

O 

1  o  o  o 

!  9  «  s? 

half  angle 
half  angle 

□ 

> 

half  angle 

x 

ILLUS 

TRATIVE  ONLY 

□ 

4 

tm* 

200,000-ft  alt 

itude ,  turbulent 

_ 1 

_ 

i 

i 

8 

0 

ILLUSTRATIVE  ONLY 

_ 

200,000-ft  altitude,  laminar 

fc 

O 

□ 

i 

8 

A 

O 

ILLUSTRATIVE  ONLY 

L-)l - 

250 ,000-f t  alti 

.tude ,  laminar 

10  20  30  40 

Vehicle  velocity,  thousands  of  feet  per  second 


50 


Figure  7.6.-  Plots  of  E^/p  versus  vehicle  velocity  for  wedges  at  zero 
angle  of  attack.  For  purposes  of  illustration  only.  (Improper  ioniza¬ 
tion  frequency  used.  See  Section  7.1.) 


Breakdown  field  strength,  volts/ cm 


b,  cm 


V  -  VQ  (1  -  e”by)  E  -  Ege  6 ^ 

O  b  ■  0  corresponds  to  infinite  boundary-layer  thickness 
or,  equivalently,  no-flow  breakdown 

0  b  ■  1.86  cm-1  is  the  actual  value  of  b  for  the 
boundary  layer  at  these  conditions 

A  b  ■  €  “  0.54  cm”1 

i 

0  b  -  lOe 

i 

Ob"  lOOe 


Figure  7.7.-  Effect  of  varying  velocity  profile  thickness 
with  velocity  outside  boundary  layer  and  all  other  param- 
eterg  held  constant.  At  conditions  of  the  flow  field  of 
a  10°  half-angle  cone  at  200,000-foot  altitude,  zero  angle 
of  attack,  and  35 ,000-f t/sec  velocity  at  a  point  5.73  feet 
downstream  of  the  vertex.  For  purposes  of  illustration 
only.  (Improper  ionization  frequency  used.  See  Section 
7.1.) 
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Figure  7.10.-  Aperture  breakdown  field  strength  versus  vehicle  velocity 
at  250,000-foot  altitude  for  wedges  and  cones  for  varying  half  angle 
at  zero  angle  of  attack.  The  field  intensity,  EB,  is  the  field 
in  the  antenna  aperture  at  which  breakdown  will  occur. 
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Figure  7.11.-  Aperture  breakdown  field  strength  divided  by  p* 
versus  vehicle  velocity  at  250,000-foot  altitude  for  wedges  and 
cones  of  varying  half  angle  at  zero  angle  of  attack.  The  field 
strength  Eg  is  the  field  in  the  antenna  aperture  at  which 
breakdown  occurs.  The  expression  p*  ■  (p/pQ)Z(760)  is  a  measure 
of  the  particle  density  in  the  boundary  layer. 


8.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

8.1  Summary  of  Work  Accomplished 

(1)  A  numerical  technique  was  formulated  for  obtaining  a 
solution  for  the  antenna  breakdown  equation  giving  the  field 
strength  at  breakdown  for  any  generalized  shape  of  electric  field 
or  flow-field  condition.  We  have  compared  this  formulation  with 
existing  experimental  data  for  nonflow  conditions. 

(2)  Breakdown  conditions  were  computed  for  a  typical  slot 
antenna  mounted  on  either  a  wedge-  or  cone-shaped  vehicle  under 
typical  high-velocity  flight  conditions. 

(3)  A  critical  examination  was  performed  of  the  electron 
conservation  equation  and-. the  concept  of  breakdown  for  c-w, 
multiple-pulse  and  single-pulse  conditions. 

(4)  Estimates  were  derived  for  the  material  properties  of 
the  gases  flowing  over  the  antennas  for  the  range  of  environmental 
conditions  expected  in  typical  flight  situations. 

(5)  The  inviscid  and  viscous  flow  fields  about  the  antenna 
were  calculated  for  several  vehicle  geometries  of  interest  and 
over  a  range  of  vehicle  velocities  and  altitudes.  These  results 
were  used  in  the  previously  discussed  calculations. 

8.2  Conclusions 

(1)  The  effect  of  the  vehicle  velocity  and  associated  shock 
wave  and  boundary  layer  on  the  breakdown  field  strength  is  not 
straightforward.  It  may  increase  or  decrease  the  breakdown  field 
as  follows: 

(A)  The  vehicle  velocity  tends  to  increase  the  breakdown 
field  strength  over  static  breakdown  by  (a)  increas¬ 
ing  the  density  of  the  gas  near  the  antenna  and 
(b)  sweeping  out  electrons  from  the  critical  region 
near  the  antenna.  The  second  effect  is  usually 
smaller  than  the  first. 
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(B)  The  vehicle  velocity  tends  to  decrease  the  break¬ 

down  field  strength  over  static  breakdown  by 
inducing  ionization  in  the  gases  above  the  antenna, 
changing  the  diffusion  from  free  to  ambipolar  and 
possibly  increasing  the  conductivity  of  the  gas  to 
the  point  where  appreciable  absorption  of  the 
energy  in  the  electric  field  results. 

(C)  The  vehicle  velocity  may  either  increase  or  decrease 

the  breakdown  field  strength  due  to  the  induced 
presence  of  ablation  products  or  the  formation  of 
new  air  species,  each  of  which  may  have  different 
ionization  and  attachment  rates  than  the  original 
air  species  and  which  may  alter  the  effective  dif¬ 
fusion  and  ionization  coefficients  of  the  air. 

In  the  presence  of  a  low  ionization  density  such  that  ambipolar 
diffusion  exists  already,  increasing  the  vehicle  velocity  usually 
tends  to  increase  the  breakdown  field  strength. 

Before  the  relative  merits  of  placing  the  antenna  in  the  high- 
pressure  upstream  regions  or  downstream  on  the  body  can  be  properly 
assessed,  one  must  take  note  of  the  possible  downstream  effect  of 
turbulence  on  electron  diffusion.  If  turbulence  increases  the 
electron  diffusion  sufficiently,  a  downstream  turbulent  position 
may  well  be  best  for  the  antenna.  We  see  that  gas  density  and 
pointing  of  the  antenna  pattern  favor  the  forward  position,  whereas 
turbulent  diffusion,  high  flow  velocity,  and  lower  background  ion¬ 
ization  favor  the  downstream  turbulent  position.  Further  extensive 
calculations  using  the  methods  developed  in  this  report  would  be 
necessary  to  assess  relative  merits  at  various  altitudes  and  veloc¬ 
ities,  or  for  a  given  trajectory.13 


13Note  that  this  discussion  concerns  only  antenna  breakdown.  The 
weak-field  plasma  interaction  usually  considered  may  introduce 
other  considerations  for  antenna  placement. 
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(2)  The  near-field  zone  of  the  antenna  is  the  important 
region  in  determining  breakdown  and  the  near  zone  field  strength 
decreases  rapidly  with  distance  from  the  antenna.  Therefore,  any 
scheme  by  which  the  effective  field  strength  at  the  edge  of  the 
flow  field  is  reduced,  such  as  withdrawing  an  antenna  below  the 
vehicle  surface  and  using  a  large  dielectric  cover,  will  greatly 
increase  the  breakdown  field  strength,  allowing  large  increases 
in  transmitted  power  without  breakdown. 

(3)  As  pointed  out  in  Section  2  the  breakdown  field  strength 
level  may  be  increased  by  using  pulses  rather  than  c-w  and  reducing 
the  pulse  widths  as  much  as  possible.  The  use  of  pulsed  systems 

is  thus  to  be  desired  if  other  considerations  permit. 

(4)  The  necessary  macroscopic  data  for  hot  air  and  for  hot 
air  with  impurities  is  almost  completely  unknown  and  the  absolute 
reliability  of  estimates  of  these  data  is  unknown  although  the 
data  used  in  this  report  are  believed  to  be  reasonable. 

8.3  Recommendations 

(1)  Additional  cases  should  be  computed  for  the  conditions 
considered  in  this  report  as  well  as  others,  using  more  terms 
for  the  expansions  in  the  z-direction. 

(2)  Controlled  experiments  similar  to  those  suggested  in 
Reference  12  should  be  performed  in  order  to  better  fix  the  macro¬ 
scopic  coefficients  necessary  to  compute  the  breakdown  phenomenon. 

(3)  Controlled  experiments  on  reentry  vehicles  should  be 
performed  in  order  to  measure  breakdown  field  strengths  (power 
levels)  as  a  function  of  vehicle  shape,  altitude,  velocity,  r.f. 
frequency,  etc.  Data  available  from  past  flights  should  be  col¬ 
lected  and  compared  with  computations  performed  using  the  techniques 
described  in  this  report. 
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APPENDIX  A 


THE  ELECTRON  KINETIC  EQUATION 
UNDER  FLOW  CONDITIONS 

In  Reference  12,  the  kinetic  equation  for  electrons,  as 
derived  in  References  1,  13,  14,  17,  18,  or  19,  is  generalized  for 
the  case  of  a  coordinate  system  in  which  the  gas  is  flowing.  It 
is  found  that  one  effect  of  the  flow  is  to  drag  the  electrons 
along  with  the  flow,  giving  rise  to  extra  collision  terms.  How¬ 
ever,  it  is  pointed  out  in  Reference  12  that  one  can  avoid  the 
consideration  of  these  terms  entirely  by  first  adopting  a  coordi¬ 
nate  system  that  moves  with  the  flow,  one  for  which  these  extra 
terms  vanish.  The  equation  is  then  essentially  the  same  as  derived 
by  the  previously  mentioned  authors,  and  we  can  then  make  the 
transformation 

7ne«  "  ?old  +  7  (‘old  *  4o) 

Slew  ™  told  ”  fco 

where  V  -  V(y)lc  is  the  flow  velocity.  Here  we  have  not  changed  the 
electron  velocity  coordinates  so  that  the  electron  velocity  dis¬ 
tribution  still  describes  the  electron  random  velocity,  c  ■  v  - 
but  in  a  stationary  coordinate  system.  The  only  change  in  the 
equation  due  to  the  transformation  is  the  appearance  of  a  V  •  VFe 
term. 

Let  us  describe  very  briefly  the  derivation  of  the  kinetic 
equation.  It  is  assumed  that  the  velocity  dependence  of  F@>  the 
density  in  phase  space  of  electrons,  may  be  expanded  in  spherical 
harmonics  about  some,  as  yet  arbitrary,  polar  axis.  A  good  argu¬ 
ment  that  terms  higher  than  the  first  two  may  be  neglected  may  be 
found  in  Reference  17.  Then  we  include  only 

F  ■  F  +  "c  •"? 
e  o  i 


(A.l) 
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when  in  the  moving  coordinate  system.  This  then,  is  substituted 
into  the  kinetic  equation.  After  discarding  collision  terms  of 
order  (ny/M)2,  the  two  linearly  independent  equations  of  lowest 
order  become 
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where  q  gives  the  rate  at  which  electrons  appear  or  reappear  at 
a  particular  velocity  after  inelastic  collisions. 

These  equations  contain  two  more  assumptions.  First,  it  is 
evident  that  charge-charge  collisions  have  been  neglected.  This 
assumes 
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en 
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since  an  electron  can  transfer  on  the  average  only  2m/M  of  its 
energy  in  a  single  collision  with  a  heavy  particle.  This  restric¬ 
tion  may  be  lifted  and  the  analysis  extended  to 

ven  »  vea  <A'5a) 

with  the  treatment  in  Reference  17.  The  results  of  this  extension, 
which  modifies  the  equations  only  slightly,  are  contained  in 
Section  3.1.1. 

Another  assumption  contained  in  Equations  (A. 2)  and  (A. 3)  is 
that  E  is  independent  of  Fe<  This  assumption  precludes  the 
possibility  of  plasma  resonance  oscillations.  The  assumption  is 
valid  if  the  applied  frequency  satisfies 
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a)  »  0)p  -  2ir  ^  J  -  2tt  (8.97x10s) 

At  A  ■  1  cm,  this  gives 

np <<  (r^~) cm~3  (A,5b) 

This,  in  practice,  is  far  less  restrictive  than  Equation  (A. 5a)  as 
expressed  in  Equation  (3.3).  However,  it  does  cause  some  diffi¬ 
culty  in  treating  the  single  pulse  breakdown  criteria  (see  Section 
2.3.2)  . 

Let  E  ■  Eq  cos  out  and  let  Fq  be  expanded  in  a  Fourier 
series 
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where  the  coefficients  may  be  functions  of  time  (i.e.,  relaxing), 
but  (dFos/dt)/FoS  «  a/2i r.  Substitution  of  this  expansion  into 
Equation  (A. 3)  and  then  solving  the  resulting  equation,  gives 
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Here  we  have  neglected  transients  in  the  solution  for  ? ,  which 

behave  as  exp(-  v  t) .  This  assumes  then,  that  (3f  s/dt)/F  8  «  v 

en  o  o  en 
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(if  it  were  not,  one  could  neglect  the  time  dependence  entirely  as 

relaxation  would  occur  quickly)  and  that  1/v  i-8  a  negligible 

?en 

■  0  at  t  ■  0. 

Substitution  of  Equations  (A. 7)  ,  (A. 6)  ,  and  E  ■  Eq  cos  cot 
into  Equation  (A. 2) ,  and  separation  of  linearly  independent  terms 
gives  for  the  equation  governing  Fq° 


(A. 8) 


The  terms  involving  F  1  are  generally  very  small  compared  to  the 
similar  terms  in  F  °  and  are  usually  neglected,1  These  terms 
are  due  to  pushing  the  alternating  component  of  the  electron  density 
back  and  forth  and  the  diffusion  away  from  this  alternating  component. 
This  alternating  component  should  be  negligible  compared  to  the 
average  electron  density,  except  near  plasma  resonance  (i.e.,  where 
Eq.  (A. 5b)  is  not  satisfied) . 

The  term  in  F  2  is  a  heating  term  due  to  second  harmonic 
oscillations.  Phenomenological  attempts  to  estimate  the  effect  of 
this  term  seem  to  indicate  that  it  is  equivalent,  at  most,  to 
increasing  the  Eeff/P  by  6  volts/ cm-mm  Hg.  This  largest  correction 
occurs  at  pX  >  5000  cm-mm  Hg  and  is  negligible  for  pX  ~  1  cm-mm  Hg 
or  less  (see,  e.g.,  Section  3.2.4). 

1  Regarding  the  neglect  of  F  1 , 
of  page  407.  ° 


see  also  Reference  13,  bottom  half 
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With  neglect  of  these  terms  and  the  subsequent  transformation 
of  the  position  and  time  variables ,  as  mentioned  above ,  Equa¬ 
tion  (A. 8)  becomes  Equation  (3.3).  Also,  substitution  of  Equa¬ 
tion  (A. 7)  into  (A.l)  gives  Equation  (3.1),  if  only  first  terms 
are  retained. 

For  the  d.c.  case  under  the  same  approximations  (which 
includes  (dF^/dt)/^  «  ve>1)  »  Equation  (A. 3)  has  the  solution 
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Substitution  of  this  into  Equation  (A. 2)  gives  the  same  equation 

as  in  the  a.c.  case  with  E .  2  replacing  E  2/(l  +  oo2/v  2) 

_  dc  rms  en 
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APPENDIX  B 


A  VARIATIONAL  PRINCIPLE  FOR  EIGENVALUES 


In  Reference  26,  page  1109,  it  is  pointed  out  that  if  r,  Cl, 
f,  and  ft  are  linear  differential  operators  such  that  they 
satisfy 
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Then  the  condition  that,  for  a  definite  pair  of  functions  0  and  0 , 
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does  not  change  with  first-order  changes  in  0  and  0  demands 
that  0  satisfy 

T0  -  k2%  (B.4) 


and  0  satisfy 

T0  -  k2fi0 
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so  long  as  the  functions  0  and  0  (and  thus  the  first-order 
changes  in  0  and  0  as  well)  satisfy  Equations  (B.l)  and  (B.2). 
Such  a  function  0  is  generally  written  0  and  called  the  adjoint 
of  0.  Also,  P  is  the  adjoint  operator  of  r  and  Cl  is  the 
adjoint  operator  of  H. 

Let  us  apply  these  facts  to  our  particular  problem.  Consider 
Equation  (5.1) 


3-2 


-  k2  « 


J  %  V2  -  V  •  (if  f  dx 


/ 


(B.6) 


fGrp  dT 


Let  us  take  the  first-order  variation 

6(k2)  =  k2  (^  +  bip,  ip  +  bip)  -  k2(^,  ip) 

where  rp ,  f,  bip ,  and  6^  all  vanish  on  the  boundaries  of  the 
region , 1  then 
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Note  that 
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where  vector  identities  and  Gauss'  integral  theorem  have  been  used. 
This  surface  integral  vanishes  because  ip  and  6  both  vanish  on 
the  boundaries.  Then  by  substituting  Equation  (B.8)  into  (B.7), 
using  (B.6)  again,  and  setting  6 (k2)  equal  to  zero,  one  obtains 


xThe  boundaries  here  are  y  ■  0,  y  ■  »,  z  ■  +_  «,  x  ■  +_  oo. 
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If  this  expression  is  to  vanish  for  all  possible  changes  in 
f  and  the  two  integrands  must  vanish.  This  gives  Equations 

(2.9)  and  (5.2). 
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DEVELOPMENT  OP  k2  WITH  ORTHOGONAL  POLYNOMIALS 


For  development  of  II 


some  identities, 
we  have 


From  Reference  27, 
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In  addition,  one  can  derive  from  Equations  (C.la)  and  (C.lb) 
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For  the  Laguerre  polynomials  we  will  need  the  identities  176.4a  and 
b  of  Reference  27,  namely,1 
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here  is  Ln(y,  a,  3)  in  Reference  27. 
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In  addition  one  can  derive2  from  Equation  (5.10a) 
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Using  Equations  (C.lb)  and  (C.lc)  there  follows: 
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Similarly  by  use  of  Equations  (C.2b)  and  (C.2c)  there  follows: 
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From  Equations  (C.3),  (C.^),  and  (5.9)  there  follows  immediately: 
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from  these  two  results  Equation  (C.2c)  follows. 
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where  QjJ||nn,  is  defined  by  Equations  (5.12),  (5.15),  and  (5.16) 
Thus  with  Equation  (5.8) 
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For  evaluation  of 
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use  of  Equation  (C.la)  and  (C.lb)  gives  immediately 
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Similarly,  it  immediately  follows  (without  any  of  the  identities 
of  Eqs.  (C.l)  or  (C.2))  that 
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APPENDIX  D 


REDUCTION  RELATIONS  FOR  THE  Iy,T) 

m,n 


By  use  of  Equation  (C.lb)  we  have 
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or  integrating  the  last  term  by  parts  and  using  Equation  (C.la) 
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From  this,  form  the  following  three  expressions:  one  by  taking 
m  ■  n  +  2;  another  by  taking  m  ■  n  +  2  and  then  increasing  n 
by  one  (n  -*  n  +  1) ;  a  third  by  taking  m  ■  n  and  then  increasing 
n  by  one.  From  these  three  expressions  eliminate  1^,’ ^  ,  and 

v  „  n+3  ,n+i 

In+2 ,n  obtaining 
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If  now  n  is  decreased  by  three  (n  -*■  n  -  3)  one  obtains  Equa¬ 
tion  (5.17).  To  obtain  Equation  (5.18)  take  m  ■  n  +  2v  in 
Equation  (D.l). 
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By  use  of  Equation  (C.2a)  we  have 
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To  obtain  an  expression  with  which  we  can  eliminate  J  from 
Equation  (E.l),  we  use  Equation  (C.2c)  twice  -  first  with  yLa*2 
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and  then  instead  with  yLa’2  -  add  the  results  and  divide  by  two 
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then  solving  for  J 
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Substitution  from  Equation  (E.2)  for  J®**]  into  Equation  (E.l) 
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gives 

-  t|  J°  ’  ^  “a  J**?  ^  -a(m+n+4-c)Ja,T^ 

1  m,n  nrt-i  ,n-i  m,n-i 


+  (2n  +  l)(a  -  T))  j“*2  +  (n2  -  1)  (a  -  n)  (E.3) 
1  m,n-i  •  m,n-2 

Examination  shows  that  this  reduction  relation  is  still  not 
adequate.  It  can  be  further  developed  by  use  of  Equation  (C.2a) . 
Note  that 


l£’2  L“’2  -  La’2 


(-  ay  +  2n  +  1)  l“’2  +  (n  -  1)  (n  +  1)  l“,2J 


"  L^,Ln  ,  +  (2n  +  1)  LL  -  (n2  -  1)  LL 
nvr  i  n— i  in  n— 1  in  n— 2 

-  (2m  +  3)  LL  +  m(m  +  2)  L  L 

m  n-i  m— i  n— 1 

Then  solving  for  Lmf1Ln_1»  multiplying  by  e“(a+Tl)y  yc  and 
integrating,  one  finds  from  this  last  expression 

n  i  "  +  2(m  -  n  +  1)  +  (n2  -  1) 

mri , n— i  m,n  m,n— i  m,n— z 


-  m(m  +  2)  , 

m-i } n*i 


(E  .4) 


Then  use  this  to  eliminate  J  .  from  Equation  (C.3)  and 

iTH-i  ,n-i 

gather  like  terms.  This  results  in  Equation  (5.19). 
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